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Executive  Summary 

This  report  details  the  development  of  mathematical  models  that  describe  the  effects  of 
acute  radiation  exposure  on  hematopoiesis.  The  models  were  developed  from  an  existing 
structure  and  modified  so  that  they  could  be  integrated  with  effects  on  hematopoiesis  from 
other  injuries,  such  as  thermal  burns.  Experimental  data  provided  insight  on  many  model 
parameters;  other  parameters  were  fit  to  case  study  data.  The  final  models  were  validated 
against  additional  case  study  data.  Challenges  were  encountered  in  using  case  study  data 
for  model  optimization  and  validation  due  to  individual  variability,  uncertainty  in  radia¬ 
tion  dose  estimates  in  human  accident  data,  and  potential  impact  of  treatment.  However, 
models  were  developed  that  capture  the  major  trends  in  circulating  platelets,  granulocytes, 
and  lymphocytes  after  radiation  exposure.  The  thrombocyte  model  demonstrates  a  dose- 
dependent  decline  in  platelets  and  subsequent  recovery.  Likewise,  the  granulocyte  model 
illustrates  an  initial  decline  in  granulocyte  count,  a  dose-dependent  abortive  rise,  and  sub¬ 
sequent  recovery.  Lymphocyte  model  predicts  dose-dependent  lymphocyte  decline  and  slow 
recovery.  The  outcome  of  our  work  is  models  that  describe  the  dose  response  of  circulating 
platelets,  granulocytes,  and  lymphocytes  as  a  function  of  time  after  irradiation.  The  infor¬ 
mation  provided  is  valuable  in  understanding  the  physiological  effects  of  radiation  exposure 
and  the  time  course  of  these  effects.  In  addition  to  combined  injury  model  development, 
future  work  could  include  the  integration  of  parameters  for  treatment,  partial  body  effects, 
and  population  variability.  These  modeling  efforts  support  more  realistic  radiological  and 
nuclear  scenario  analyses  in  large  diverse  populations,  thereby,  facilitating  preparedness  and 
medical  resource  planning. 
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Introduction 


1.1  Background 

Applied  Research  Associates,  Inc.  (ARA)  has  been  tasked  by  Defense  Threat  Reduction 
Agency  (DTRA)  to  support  their  mission  to  safeguard  against  weapons  of  mass  destruction 
(WMD).  ARA  is  supporting  this  effort  by  developing  state-of-the-art  mathematical  models 
that  predict  medical  and  performance  consequences  from  radiation  and  combined  injuries, 
thereby  enhancing  our  understanding  of  the  potential  impact  of  a  nuclear  detonation.  This 
work  will  improve  current  casualty  estimation  capabilities  through  an  interdisciplinary  ap¬ 
proach  that  integrates  experimental  data  with  mechanistic  mathematical  modeling. 

The  hematopoietic  models  presented  in  this  work  describe  hematopoiesis  after  acute  ra¬ 
diation  exposure  without  treatment.  The  main  intent  in  developing  these  models  was  to 
merge  them  with  analogous  models  of  acute  injury,  such  as  burn,  to  better  understand  the 
impact  of  combined  injury  on  pathophysiology  and  thereby  improve  mortality  estimates  of 
combined  injury  by  accounting  for  mechanistic  interactions.  Hematopoietic  model  param¬ 
eters  for  thermal  injury  have  also  been  developed;  however,  details  on  the  models  for  burn 
and  combined  injury  will  be  presented  subsequently.  In  addition  to  the  insight  gained  from 
combined  injury  modeling,  the  models  of  hematopoiesis  and  radiation  alone  provide  clini¬ 
cally  relevant  outputs  that  afford  insight  on  the  time-dependent  evolution  of  casualties  and 
hence  medical  resource  requirements  as  a  function  of  time. 

The  hematopoietic  models  developed  in  our  current  work  will  complement  the  Radiation- 
Induced  Performance  Decrement  (RIPD)  model,  which  is  a  physiologically-based  mechanistic 
model  of  acute  radiation  sickness  (ARS)  developed  in  the  1990s  (Matheson  et  al.  1998).  RIPD 
was  integrated  into  a  software  tool  used  for  planning  and  scenario  predictions  for  nuclear 
detonations,  has  since  been  updated,  and  will  be  included  in  our  broader  software  platform 
Health  Effects  from  Nuclear  and  Radiation  Environments  (HENRE).  HENRE  will  house  all 
new  models  developed  for  prompt  injury,  combined  injury,  and  other  relevant  models  as  they 
are  developed. 

1.2  Purpose 

Physiologically-based  mechanistic  models  of  the  human  response  to  irradiation  can  yield 
valuable  information  for  use  in  nuclear  disaster  preparedness  planning.  Understanding  how 
biological  systems  change  after  radiation  exposure  provides  insight  on  the  pathophysiolog¬ 
ical  processes  leading  to  mortality  as  well  as  the  evolution  of  these  processes  over  time. 
Models  that  describe  these  processes  can  be  used  to  correlate  biological  endpoints  with 
time-dependent  clinical  manifestations  and  probabilities  of  mortality.  Information  on  clin¬ 
ical  endpoints  can  help  medical  planners  determine  when  and  under  what  circumstances 
medical  resources  might  be  needed.  For  example,  the  time  at  which  platelet  counts  decrease 
below  a  critical  threshold  can  suggest  when  platelet  transfusions  will  be  needed,  and  when 
analyzed  on  a  population  level,  how  many  transfusions  for  a  particular  scenario  might  be 
required. 

This  report  presents  our  work  on  models  describing  hematopoiesis  after  acute  radiation 
exposure.  Hematopoiesis  is  the  process  by  which  blood  cells  are  generated  from  hematopoi- 
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etic  stem  cells  (HSCs)  in  the  bone  marrow.  The  cellular  components  of  blood  play  essential 
roles  in  nourishing,  protecting,  and  healing  the  human  body.  Due  to  their  fast  prolifer¬ 
ation  rate,  bone  marrow  hematopoietic  progenitors  are  highly  radiosensitive.  Thus,  the 
hematopoietic  system  shows  a  marked  response  to  irradiation.  The  hematopoietic  syndrome 
of  ARS  has  been  well-documented  and  characterized  (Flicdner,  Friesecke,  et  al.  2001).  The 
hematopoietic  response  after  radiation  has  been  used  to  estimate  radiation  dose  after  acci¬ 
dental  over-exposures,  develop  prognoses,  and  guide  treatment  decisions. 

Our  work  on  physiologically-based  mechanistic  models  describing  the  effects  of  radiation 
exposure  in  humans  on  thrombopoiesis,  granulopoiesis,  and  lymphopoiesis  is  detailed  in  this 
report.  Thrombopoiesis,  granulopoiesis,  and  lymphopoiesis  correspond  to  the  processes  of 
formation  of  platelets  (also  known  as  thrombocytes),  granulocytes,  and  lymphocytes,  respec¬ 
tively.  Granulocytes  and  lymphocytes  are  both  white  blood  cells  which  are  highly  involved 
in  the  immune  response.  Platelets,  which  have  a  less  well  understood  role  in  immunity, 
are  also  involved  in  hemostasis  and  help  with  coagulation,  or  blood  clot  formation.  These 
models  describe  hematopoietic  cell  sensitivity  to  irradiation,  i.e.  radiation  dose-dependent 
rate  of  cell  damage  and  death,  and  the  impact  of  the  perturbed  state  on  cell  proliferation 
and  differentiation  rates.  The  methodologies  used  for  model  development,  parameterization, 
computational  algorithms,  optimization,  and  validation  are  described  in  this  report. 
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2  Methods 


2.1  Existing  Mathematical  Models 

A  number  of  mathematical  models  of  hematopoiesis  have  been  developed  previously  for  un¬ 
derstanding  blood  cell  dynamics  in  general  and  following  different  insults  and/or  treatment. 
These  include  models  of  radiation  exposure  (Grafile  2000;  Akushevich  et  ah  2011),  granulo¬ 
cyte  colony  stimulating  factor  treatment  (Vainstein  et  ah  2005),  and  chemotherapy  (Scholz 
et  ah  2010).  Russian  scientist,  Olga  A.  Smirnova,  developed  a  series  of  mathematical  models 
that  describe  cellular  dynamics  after  acute  and  protracted  radiation  exposures  (Smirnova 
2011).  These  models  are  based  on  several  decades  of  radiobiological  experiments  and  human 
case  study  data.  Dr.  Smirnova’s  series  of  hematopoietic  models  include  murine  throrn- 
bopoiesis,  granulopoiesis,  and  lymphopoiesis  and  human  thrombopoiesis  and  granulopoiesis 
(Smirnova  2011;  Smirnova  2012).  NASA  has  also  developed  radiation  models  based  off  of 
Smirnova’s  work  including  models  for  lymphopoiesis  in  humans  (Hu,  Smirnova,  et  ah  2012) 
and  granulopoiesis  in  rodents,  canines,  primates,  and  humans  (Hu  and  Cucinotta  2011). 
Since  these  models  describe  the  impact  of  radiation  on  blood  cell  kinetics,  they  served  as 
excellent  starting  points  for  our  model  development.  However,  data  and  details  on  how  some 
parameters  were  derived  were  lacking  or  embedded  in  Russian  literature.  An  attempt  was 
made  to  reconstruct  as  much  of  Smirnova’s  original  work  as  possible,  including  obtaining 
and  translating  Russian  literature  that  was  accessible.  Nevertheless,  many  details  remained 
unclear.  Because  Smirnova’s  models  were  designed  for  simulating  both  acute  and  protracted 
radiation  effects,  some  parameter  values  were  likely  determined  so  that  protracted  effects 
could  be  accurately  predicted.  The  current  work  is  focused  on  predicting  the  hematopoietic 
response  to  acute  (single,  high  dose  rate)  radiation  exposures;  we  have  not  endeavored  to 
address  protracted  doses  at  this  point. 

2.2  Model  Development 

In  this  work  parameter  values  were  tuned  to  match  acute  effects  as  accurately  as  possible. 
The  human  model  development  was  guided  by  work  done  on  establishing  non-human  primate 
(NHP)  hematopoietic  models  for  a  separate  project  (Wentz  et  ah  2014).  For  the  NHP 
model  development,  a  significant  amount  of  hematological  data  from  experimental  work 
with  rhesus  monkeys  was  made  available  to  us  from  a  clinical  partner.  This  data  was 
supplemented  with  published  data  on  rhesus  hematopoietic  kinetics  after  different  irradiation 
doses.  Key  features  of  the  hematopoietic  response  and  required  structure  of  the  models  could 
be  delineated  from  work  with  the  NHP  data  since  these  studies  involved  precise  experimental 
conditions.  Modern  literature  was  reviewed  and  experimental  data  collected  for  as  many 
of  the  hematopoietic  model  parameters  as  possible.  Based  on  the  NHP  work,  availability 
of  experimental  data,  and  published  parameter  values  (for  example,  compartment  sizes  or 
transition  rates),  the  model  structure  was  altered  as  needed  and  new  model  parameters  were 
developed. 

To  develop  human  model  parameters,  optimize  their  values,  and  validate  model  outputs, 
clinical  data  on  human  blood  cell  kinetics  after  radiation  exposure  was  required.  Therefore, 
published  data  on  radiation  accidents  and  radiation  therapy  was  collected.  Details  on  the 
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case  studies  reviewed  and  selected  are  described  in  more  detail  in  Appendix  A.  Briefly, 
data  available  on  acute,  whole  body  gamma  exposures  that  had  little  or  no  treatment  was 
prioritized.  Since  very  little  data  is  available  for  strictly  “acute,  whole  body  gamma  exposures 
with  no  treatment”,  a  number  of  other  studies  were  used  with  careful  consideration  placed 
on  assigning  relevant  doses  in  cases  involving  different  quality  radiations,  inhomogeneous 
exposures,  protracted  exposures,  and  treatment. 

2.3  Parameter  Estimation  and  Confidence  Interval  Calculations 

The  statistical  computing  environment  R  (v3.0.0,  R  Core  Team  2013)  was  used  to  generate 
numerical  solutions  of  ordinary  differential  equations  (ODEs),  estimate  model  parameters, 
and  provide  confidence  intervals  for  parameters.  Specifically,  ODEs  were  solved  using  the 
“ode”  function  from  the  deSolve  package  (Soetaert,  Petzoldt,  and  Setzer  2010),  and  param¬ 
eter  estimation  was  performed  using  the  modCost  and  modFit  tools  from  the  FME  package 
(Soetaert  and  Petzoldt  2010).  ModFit  performs  constrained  fitting  of  a  model  to  data  using 
input  including  (1)  a  function  to  be  minimized,  (2)  initial  parameter  values,  (3)  the  lower 
and  upper  bounds  to  be  used  for  each  parameter,  and  (4)  the  optimization  method. 

Parameter  estimation  can  be  formulated  as  a  nonlinear  least  squares  optimization  prob¬ 
lem,  and  the  structure  of  the  least  squares  objective  function  can  be  leveraged  to  create  alter¬ 
natives  to  Newton’s  method,  such  as  the  Levenberg-Marquardt  Method  (More  1978;  Heath 
1997).  The  Levenberg-Marquardt  Method  can  be  viewed  as  a  combination  of  the  Gauss- 
Newton  Method  and  the  Steepest  Descent  Method,  and  it  can  provide  a  robust  and  efficient 
way  to  locate  optimum  parameter  values.  The  modFit  function’s  Levenberg-Marquardt 
option  was  thus  used  as  part  of  our  parameter  estimation  methodology. 

The  objective  function  /  we  provided  to  modFit()  maps  a  given  parameter  vector  to 
a  residual.  This  mapping  includes  the  following  steps: 

1.  The  ODE  simulation  is  performed  at  each  radiation  dose  level. 

2.  The  output  of  each  ODE  simulation  and  the  observed  data  is  used  as  input  for  modCost 
which  calculates  and  returns  the  residuals. 

3.  The  residuals  from  each  simulation  are  combined  across  doses  and  returned. 

modFit  then  uses  /(£*,)  and  the  Levenberg-Marquardt  Method  to  generate  the  next  ap¬ 
proximate  parameter  vector  C,k+i-  These  steps  are  then  repeated  to  generate  additional 
approximate  solutions  C.k+2,  £fc+3,  •  •  •  until  a  convergence  test  is  passed.  Finally,  we  note 
that  the  optimizations  were  performed  on  normalized  data  in  which  each  data  value  was  di¬ 
vided  by  a  baseline  value.  This  allows  for  comparison  of  residuals  between  datasets  without 
inadvertently  weighting  one  dataset  more  than  another. 

To  calculate  parameter  confidence  intervals,  the  Markov  Chain  Monte  Carlo  (MCMC) 
method  with  the  adaptive  Metropolis  algorithm  was  used  (Laine  2008;  Haario  et  al.  2001). 
Given  an  initial  parameter  vector  £j,  an  estimate  for  the  parameter  covariance  matrix,  and 
the  same  function  /  as  described  above,  the  parameter  space  was  randomly  sampled  to 
generate  a  new  parameter  vector  £j+i.  The  likelihood  L  of  the  new  parameter  vector  £j+ 1 
was  determined,  and  if  L(£j+ 1)  >  L(£j),  £J+i  was  chosen.  If  L(£j+ 1)  <  L(£j),  was  selected 
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with  probability  a,  where  a  =  P(£j+i)/P(£j)  and  P  is  the  probability  of  a  given  parameter 
vector.  These  calculations  were  performed  in  R,  using  modMCMC  which  is  also  a  function 
in  the  FME  library.  To  determine  whether  the  models  were  identifiable  the  collincarity  value 
was  determined  using  the  collin  function  from  the  FME  package  in  R  (see  Brun  et  al.  2001). 


7 


This  page  intentionally  left  blank. 


3  Overview  of  Hematopoietic  Models 


3.1  Model  Structure 


To  model  the  complex  process  of  hematopoiesis,  several  simplifications  that  preserve  key 
system  characteristics  are  adopted.  Cells  are  grouped  into  compartments  based  on  degree  of 
maturity,  differentiation,  and  location  in  the  body  (ex.  bone  marrow,  tissue,  or  circulation). 
The  primary  cell  compartments  are  mitotic  bone  marrow  precursor  cells  (Ad),  post-mitotic 
bone  marrow  precursor  cells  (X2),  and  the  fully  differentiated  cells  (granulocytes,  lympho¬ 
cytes,  or  platelets)  in  the  blood  (A3).  The  granulopoiesis  model  includes  an  additional 
compartment  (Ad)  that  represents  granulocytes  in  the  tissue. 

Ordinary  differential  equations  are  used  to  describe  the  transition  between  the  develop¬ 
mental  stages: 


(1) 

(2) 

(3) 


X\  =  Bx  1  —  Cx  1 
x2  =  Cx  1  —  Fx2 
x3  =  Fx2  -  Ex3 


where  X{  represents  the  change  in  cell  i  (i  =  1, 2, 3)  concentration  with  time.  X\  cells 
reproduce  at  rate  H,  which  is  dependent  on  the  concentration  of  the  mediator  in  the  system 
(see  Section  3.2).  C  and  F  are  rates  of  transition  between  states.  These  rates  can  be 
dependent  on  the  state  variables,  and  thus  describe  feedback  mechanisms.  E  is  the  rate  of 
normal  x3  cell  death  or  clearance  from  the  circulation. 

3.2  Stimulatory  Hematopoiesis  Mediator 

All  the  models  include  a  generic  mediator  /  as  part  of  a  negative  feedback  system  to  regulate 
the  production  of  HSCs.  In  Smirnova’s  models,  this  mediator  is  considered  to  be  inhibitory. 
That  is,  increased  mediator  concentration  directly  reduces  the  X\  repopulation  rate,  and  as 
the  number  of  cells  in  the  model  increases,  more  of  the  mediator  is  produced.  However,  we 
have  changed  the  biological  interpretation  of  the  mediator  from  inhibitory  to  stimulatory 
in  order  to  relate  the  mediator  to  the  action  of  known  hematopoietic  stimulators  such  as 
thrombopoietin  (TPO)  and  granulocyte  colony-stimulating  factor  (G-CSF).  Thus,  in  our 
version  of  the  model,  the  negative  feedback  mechanism  is  less  direct:  as  the  number  of 
hematopoietic  cells  increases,  the  stimulatory  mediator  decays  more  rapidly,  leading  to  less 
stimulation  and  a  reduced  x\  repopulation  rate.  However,  in  both  versions  of  the  model 
(inhibitory  mediator,  stimulatory  mediator),  the  mediator  is  assumed  to  be  in  a  quasi¬ 
steady  state  (i.e.,  dl/dt  ~  0);  therefore,  as  explained  below,  for  this  particular  study  our 
change  of  perspective  does  not  affect  the  model’s  predictions. 

I11  the  original  model,  the  mediator  /  is  described  as  an  inhibitor  of  the  X\  repopulation 
rate.  The  concentration  of  /  increases  as  the  concentration  of  each  cell  population  increases 
according  to  the  following  equation: 


—  —  Gi(x  1  +  92x2  +  #3X3)  —  HJi 


(4) 
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Here  the  subscript  i  refers  to  parameters  in  the  original  model  that  describe  the  inhibitory 
mediator  (the  subscript  s  will  be  used  to  denote  parameters  that  describe  the  stimulatory 
representation).  Equation  4  can  be  considered  “fast”  relative  to  the  change  in  cell  concen¬ 
trations  with  time  (iq,  i;2,  £3)  and  thus  the  approximation  =  (Gi/ H{)(xi  +  d2X2  +  #3X3)  is 
used.  The  rate  of  aq  repopulation,  Bi:  is  given  by  the  following  equation: 


1  +  h/Ki 


a 

1  +  A  (2 q  +  O2X2  +  83X3) 


(5) 


where  A  =  and  a  ^ie  maximum  rate  of  repopulation. 

To  change  the  description  of  the  mediator  from  inhibitory  to  stimulatory,  the  mediator’s 
relation  to  the  concentration  of  cells  in  the  system  is  altered.  Rather  than  having  cell 
concentrations  affect  the  rate  of  mediator  generation,  cell  concentrations  affect  the  rate  of 
mediator  decay.  The  following  equation  describes  this  new  system: 

— —  =  Gs  —  Hs(x  1  +  92X2  +  @3X3)  Is  (6) 

at 


Here  the  mediator  is  produced  at  a  constant  rate  Gs.  It  is  then  removed  from  the  system  at 
a  rate  that  is  dependent  on  cell  concentration. 

The  rate  of  aq  repopulation  is  changed  so  that  K  is  a  stimulatory  rate  constant  rather 
than  an  inhibitory  rate  constant  and  is  given  by  the  following  equation. 

a  a  a 

s  ~  1  +  Ks/Is  ~  l  _|_  K‘Hs(x i+e2x2+e3x3)  -  l  +  /3s(aq  + 03X2  +  03X3) 

Cr  s 

where  A  =  757"  • 

The  equations  for  Bi  and  Bs  are  the  same  except  for  the  constants  A  and  A-  However, 
after  scaling  the  models  using  their  equilibrium  points,  the  models  are  the  same  (i.e. ,  A  and 
A  do  not  appear  in  the  nondimensional  form  of  the  models).  Thus,  using  a  quasi-steady- 
state  approximation  for  /,  the  mathematical  models  are  identical,  but  the  definition  of  the 
mediator  action  has  changed.  However,  restructuring  the  model  based  on  a  stimulatory 
mediator  affords  a  mechanism  for  integrating  treatment  with  hematopoietic  stimulators, 
such  as  G-CSF,  into  future  radiation  effects  models. 


3.3  Cell  Survival  Curves  and  Damaged  Compartments 

The  radiosensitivity  of  a  given  compartment  Xt  depends  on  the  cell  line  being  modeled.  I11 
the  thrombopoiesis  model,  cells  in  the  X\  and  X2  groups  are  sensitive  to  radiation;  in  the 
granulopoiesis  and  lymphopoiesis  models,  all  cell  groups  (Xj,  1  <  %  <  4  for  granulopoiesis,  1 
<  i  <  3  for  lymphopoiesis)  are  sensitive  to  radiation. 

If  cell  group  Xt  is  radiosensitive,  cells  in  the  group  can  become  damaged  {Xf  compart¬ 
ment).  Damaged  cells  die  within  days  at  rate  /a.  In  addition,  for  granulopoiesis,  cells  in 
the  Ad  group  can  become  weakly  damaged  (X™d  compartment).  The  weakly  damaged  state 
provides  a  way  to  model  the  abortive  rise  phenomena,  observed  5  to  30  days  after  radiation 
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exposure,  in  which  a  second  transient  increase  in  granulocyte  cell  count  occurs  (Bond  et  al. 
1965).  Weakly  damaged  cells  reproduce  and  mature  like  undamaged  cells  for  a  specified  time 
interval  A tae  after  exposure.  After  this  time  interval,  the  weakly  damaged  cells  and  their 
progeny  die  off  at  rate  rj. 

In  Smirnova’s  original  model,  cell  damage  is  modeled  using  the  one-target-one-hit  theory 
of  cell  damage  in  which  the  specific  rate  of  damage  is  proportional  to  the  radiation  dose 
(Joiner  2009).  After  acute  radiation  of  dose  D,  the  proportion  of  undamaged  or  surviving 
cells  is 


S  =  e~DiD°  (8) 

However,  this  model  is  limited  in  its  ability  to  capture  the  dose-response  characteristics  of 
cell  survival.  Therefore,  we  have  implemented  a  more  general  cell  survival  model,  the  multi- 
target-single-hit  model  (Joiner  2009).  In  this  model,  the  surviving  fraction  of  cells  after 
acute  radiation  of  dose  D  is 


S  =  1  -  (1  -  e-D/Di)ni 


(9) 


where  7i;  is  the  average  number  of  hits  required  to  damage  a  cell  in  compartment  Xt  and  D® 
is  the  dose  that  results  in  (on  average)  one  hit  per  target.  As  expected,  when  n*  =  1  this 
model  simplifies  to  the  one-tar  get-one- hit  model.  Thus,  for  each  radiosensitive  cell  group, 
there  are  two  parameters  that  determine  the  proportion  of  cells  that  are  undamaged:  D® 
and  Hi. 

For  the  granulopoiesis  cell  lineage,  an  x\  cell  can  become  either  fully  damaged  or  weakly 
damaged.  To  simulate  this,  we  define  two  different  types  of  radiation  hits,  a  full-damage  hit 
and  a  weak-damage  hit.  In  order  for  an  x\  cell  to  become  fully  damaged,  it  must  receive 
n i  full-damage  hits,  while  all  other  cells  that  receive  ri\  hits  become  weakly  damaged.  The 
fraction  of  damaged  cells  that  are  fully  damaged  following  acute  radiation  dose  D  is  given 
by  1/(1  +  $),  where 


d 


(1  _  e -D/D^m  _  (i  _ 

(1  _  e-°/fli00)" i 


(10) 


and  D®00  is  the  dose  that  results  in  (on  average)  one  full-damage  hit  per  target.  By  definition 
D®00  >  £>°. 
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4  Thrombopoiesis  Model 

4.1  Conceptual  Model 

Thrombopoiesis  is  the  process  of  platelet  formation.  Platelets  are  responsible  for  the  coag¬ 
ulation  of  blood  and  are  a  source  of  growth  factors  which  regulate  cell  growth  and  division. 
Platelet  generation  is  initiated  by  self-renewing  HSCs  in  the  bone  marrow.  These  stem  cells 
differentiate  and  through  a  series  of  mitotic  divisions  produce  megakaryocytes  (MKs),  the 
precursors  of  platelets.  As  MKs  mature,  their  ploidy  increases  through  a  series  of  endomi- 
toses.  Once  mature,  MKs  produce  platelets  which  enter  the  blood  stream. 

Similar  to  the  baseline  model  structure,  the  thrombopoietic  model  consists  of  three  com¬ 
partments.  The  X1  compartment  represents  mitotic  precursors  in  the  bone  marrow,  ranging 
from  HSCs  to  megakaryoblasts.  The  X2  compartment  represents  MKs  in  the  bone  marrow, 
and  the  X:i  compartment  represents  platelets  in  circulation.  Each  MK  produces  thousands  of 
platelets  resulting  in  an  amplification.  The  number  of  platelets  produced  per  MK  is  denoted 
a. 

Following  radiation  exposure  in  humans,  there  is  a  delay  before  the  peripheral  platelet 
counts  start  to  decline  (Bond  et  al.  1965).  This  delay  is  likely  due  to  the  time  needed  for  the 
early  progenitors  affected  by  radiation  to  mature.  In  order  to  simulate  this  delayed  effect, 
sub-compartments  were  added  to  the  MK  compartment.  A  cell  entering  the  X2  compartment 
travels  through  each  of  the  subcompartments  before  producing  platelets.  Furthermore,  the 
first  half  of  the  MK  sub  compartments  are  labeled  as  immature  MKs  and  the  second  half 
are  labeled  mature  MKs.  The  reason  for  this  differentiation  is  for  feedback  purposes  and  is 
described  in  Section  4.1.2.  The  platelet  compartment  is  also  divided  into  subcompartments 
to  generate  more  biologically  realistic  transit  time  distributions  (Murphy  and  Francis  1971). 
Figure  4.1  illustrates  the  structure  of  the  thrombopoiesis  model. 

4.1.1  Radiation  damage 

In  our  model,  both  mitotic  progenitors  and  MKs  are  radiosensitive  while  in  Smirnova’s  model 
MKs  are  radioresistant.  MK  radiosensitivity  was  added  because  at  high  radiation  doses  (5-8 
Gy),  the  delay  in  platelet  decline  following  radiation  exposure  does  not  occur,  suggesting 
MKs  are  becoming  damaged  (Cohn  and  Milne  1956).  Thus,  we  included  MK  radiosensitivity 
in  our  human  model. 

When  cells  become  damaged,  they  enter  into  only  one  damaged  state  for  each  cell  type 
(Xf  for  damaged  mitotic  progenitors  and  A|  for  damaged  MKs).  In  Smirnova’s  model  cells 
may  become  moderately,  heavily,  or  weakly  damaged.  This  level  of  detail  may  have  been 
required  for  Smirnova’s  protracted  radiation  models;  however,  none  of  the  kinetic  blood  cell 
count  data  for  acute  exposures  supported  including  a  weakly  damaged  cell  compartment,  and 
the  proportion  of  cells  entering  the  moderately  verses  heavily  damaged  cell  compartment  did 
not  significantly  affect  model  output.  Thus,  to  simplify  the  model  and  minimize  parameters, 
the  weakly  damaged  cell  compartments  have  been  removed,  and  the  moderately  and  heavily 
damaged  cell  compartments  have  been  merged  into  one  damaged  compartment. 
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Cell  transitions 


-  ■  »  Instant  radiation  damage 


■>  Feedback  effects 


Figure  4.1:  Thrombopoiesis  model  diagram 


4.1.2  Feedback  mechanisms 

The  model  has  two  feedback  mechanisms  that  are  essential  for  maintaining  and,  if  concen¬ 
trations  are  perturbed,  returning  to  equilibrium.  The  first  mechanism  involves  a  generic 
mediator,  the  level  of  which  is  regulated  by  the  concentration  of  cells  in  the  system.  The 
mediator  stimulates  the  repopulation  rate  of  mitotic  progenitors.  Because  the  action  of  the 
mediator  is  fast  with  respect  to  the  hematopoietic  cell  kinetics,  the  rate  of  repopulation  ( B ) 
can  be  written  as  a  function  of  cell  concentration  (see  Section  3.2  for  derivation). 


1  +  (3(x±d  +  xd  +  d2(x2d  +  xd)  +  93x3) 

Due  to  the  significant  stimulatory  role  of  TPO  in  thrombopoiesis  (Kaushansky  2005),  it 
is  assumed  to  represent  a  large  portion  of  the  generic  mediator.  Thus,  known  biological 
mechanisms  involving  TPO  are  used  to  justify  the  effects  of  the  generic  mediator  in  our 
model.  TPO  concentration  is  regulated  by  platelets  and  MKs  through  receptor  interactions. 
Specifically,  receptors  on  the  membranes  of  platelets  and  MKs  bind  to  and  engulf  TPO, 
removing  it  from  the  system. 

In  the  second  feedback  mechanism,  the  concentration  of  platelets  alters  the  maturation 
rate  of  immature  MKs.  Information  on  the  impact  of  platelets  on  MK  maturation  can  be 
gleaned  from  studies  on  platelet  kinetics  following  thrombopheresis,  which  is  the  removal 
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Figure  4.2:  Feedback  effect  of  platelets  x3  on  MK  maturation  time  (t2) 


of  platelets  from  an  individual  for  donation  (Sullivan  et  al.  1977).  In  these  studies,  the 
platelet  concentration  increases  at  the  expected  rate  for  the  first  3-4  days  given  known 
platelet  production  rates.  However,  an  increase  in  the  production  rate  approximately  4  days 
after  thrombopheresis  suggests  that  the  platelet  concentration  affects  the  maturation  of  early 
MKs.  TPO  is  known  to  affect  MK  maturation  as  well  as  the  early  progenitor  repopulation 
but  does  not  affect  the  formation  of  platelets  from  the  MK  cytoplasm  (Kaushansky  2005). 
For  these  reasons,  a  negative  effect  of  platelet  concentration  was  added  on  only  immature 
MK  maturation  rather  than  on  the  entire  MK  compartment.  To  divide  the  X2  compartment 
into  immature  and  mature  MKs,  the  first  half  of  the  sub-compartments  were  classified  as 
immature  X2j  while  the  second  half  of  the  sub-compartments  were  classified  as  mature  X2m- 
Thus,  the  number  of  MK  compartments  is  limited  to  being  an  even  integer. 

To  model  this  second  feedback  mechanism,  we  assume  that  the  normal  transit  time 
through  the  entire  MK  compartment  is  r0,  and  that  this  transit  time  ranges  from  Tmin  to 
Tmax ■  We  also  assume  that  (see  Fig.  4.2) 

•  the  transit  time  t2m  through  the  mature  MK  compartment  ( X2m )  is  fixed  at  To/2,  so 
that  the  transit  time  r2j  through  the  immature  MK  compartment  {X2j)  ranges  from 
T2I  ,min  7 min  to/2  to  T2i,max  —  Tmax  to/2 

•  Tmax  >  t0  >  Tmin  >  y  ( Tmin  must  be  greater  than  y  because  y  is  the  transit  time 
through  the  mature  MK  sub-compartments) 

When  the  system  is  at  equilibrium,  we  want  t2j  =  t0/2  (so  that  t2/  +  t2M  =  t0).  Moreover, 
as  x3  decreases  we  want  r2I  — >  T2jtmin,  and  as  x3  increases  we  want  t2j  — *  r2Ijmax.  With  these 
goals  in  mind,  we  define 

T21  ~  ^ i. — yij;  — ~ — ;  — ■ — ) 

,max  T~2I  ,min  ^ 
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Xmin,  Xmax ) 

<50 


•Emin  T  (^max  Emin)  3  (x  X min )  3 

1 

TO 

1 

1~max 

1 


^ ~min 


(14) 

(15) 

(16) 
(17) 


Smin  and  5max  are  respectively  the  minimum  and  maximum  possible  net  transition  rates 
through  the  entire  MK  compartment.  Note  that  the  above  assumption  that  Tmax  >  To  > 
Emin  >  -f  is  equivalent  to  the  assumption  that  <5m;n  <  5o  <  Smax  <  2 5q. 

Given  this  model,  it  is  clear  that  when  the  system  is  at  equilibrium  (and  thus  x3  =  1), 
r2 1  =  tq/ 2,  so  that  the  total  MK  transit  time  is  t2j  +  t2m  —  To  —  Sq1.  As  the  concentration 
of  platelets  approaches  zero,  r2j  approaches  (1  —  Smax/250)/Smax  =  T2ymm.  Thus,  the  total 
MK  transit  time  approaches  r2/imjn  +  t2m  =  ijmin  —  ro/2)  +  ro/2  =  rmin.  Similarly,  as  the 
concentration  of  platelets  increases,  the  total  MK  transit  time  approaches  Tmax- 

Smirnova’s  original  model  included  a  feedback  effect  of  X3  on  ploidy  levels.  However,  we 
found  that  after  adding  the  feedback  on  maturation  rate,  an  effect  of  x3  on  ploidy  did  not 
improve  model  outputs  and  was  therefore  removed  from  the  model. 

The  sources  for  transition  rates  between  compartments  and  of  cell  decay  (7,  n,  5,  Smin , 
Smax,  m,  %jj,  fj),  radiosensitivity  parameters  (D®,  n  1,  D 2,  n2),  and  other  parameters  ( a ,  A,  a, 
02 ,  6*3)  are  provided  in  Section  4.3. 


4.2  Mathematical  Model 

The  mathematical  model  of  thrombopoiesis  is  described  by  a  set  of  initial  conditions  and 
differential  equations  that  describe  the  change  in  cell  concentrations  with  time.  A  description 
of  each  cell  compartment  is  provided  in  Table  4.1. 

At  equilibrium  cell  concentrations  in  the  Ad,  X2,  and  X3  compartment  are  equal  to 
x'ld,  x2d,  and  x3,  respectively.  Following  acute  exposure,  radiosensitive  cells  immediately 
transition  into  the  damaged  compartment.  The  following  equations  define  how  cell  concen¬ 
trations  change  following  an  acute  exposure  of  dose  D  and  determine  initial  conditions  for 
the  temporal  model. 


xf  =  xf  (l  -  (1  -  e-D/D°)ni^ 

ud  /  \ 

4?,  =  d-  (1  -  (1  -  e-D^r-) 
4,  =  (1  -  (1  -  e-^r-) 


(18) 

1,  2,  ■■■,  n/2) 

(19) 

1,  2, ...,  n/2) 

(20) 
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Table  4.1:  Cell  concentrations  in  each  compartment  of  the  thrombopoiesis  model 


Variable  name  Description 

x^d  Undamaged  cells  in  the  mitotic  compartment  X\ 

x^ji  Undamaged  immature  MKs  in  the  ith  subcompartment  of  X2j 

x^Mi  Undamaged  mature  MKs  in  the  ith  subcompartment  of  X2M 

Xz,i  Undamaged  circulating  platelets  in  the  ith  subcompartment  of  X3 

xd  Damaged  cells  in  the  mitotic  compartment  X\ 

x2  Damaged  MKs  in  the  X2  compartment 

x/d  Concentration  of  cells  in  the  X%  compartment  at  equilibrium  (i=l,2) 

x3  Concentration  of  cells  in  the  X3  compartment  at  equilibrium 


X3,i  =  ~  (i  =  1,2,  (21) 

m 

xi  =  xf  (l  -  e-D/D°V'  (22) 

4  =  xf  (l  -e-D/D2)"  (23) 

Once  radiation  exposure  occurs,  a  series  of  differential  equations  are  used  to  describe 
the  compartmental  transitions,  repopulation  and  decay  rates,  and  the  feedback  mechanisms 
shown  in  Figure  4.1.  a  is  a  constant  which  represents  the  number  of  platelets  produced 
per  MK.  The  following  equations  describe  the  model  where  Eq.  18-23  are  used  to  calculate 
initial  conditions.  x,b  represents  dxi/dt. 


x?  =  Bxf  -  'yxf  (24) 

xflx  =  7 xf  -  n5x%}A  (25) 

x2i,i  =  nSx2i,i-i  -  n$x2i,i  (*  =  2,  3, ...,  n/2)  (26) 

X2M,1  =  n$x gn/2  -  n5 0^1  (27) 

X2 M,i  =  n^0%2M,i—l  ~  n^0x2M,i  (*  =  2,  3,  ...,  n/2)  (28) 

x3,l  =  anSox2M,n/2  -  ™>i>x 3,1  (29) 

x3ti  =  m'ijjx3^1  -  rmjjx3ti  (i  =  2, 3, ...,  m)  (30) 

xd  =  -^xd  (i  —  1,  2)  (31) 

where 


1  +  /3  (d1(x^d  +  xf)  +  92(x %d  +  xf)  +  03x 3) 

n/2 

ud  _  V"''  ud  |  ud 
x2  —  x2 1,i  '  x2 M,i 

i= 1 


(32) 

(33) 
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(34) 


x3  =  ^2  X3’‘ 

i= 1 

i 


<*  =  2Z<m o.  J 

Z ( X j  Xmax)  X-mifi  -\-  (x 


^mir, 


28o 

max  X 


m  I  _ 


2Sq 


\  ^  X3  ( nr  —  nr 

l  V  ^  iXj  r; 


(35) 

(36) 


4.3  Radiation  Model  Parameter  Sources 

Thrombopoiesis  parameters  were  determined  either  using  experimental  data  or  through  op¬ 
timization.  The  optimized  parameters  include  7,  <5,  A,  62,  and  D®.  The  following  gives 
a  detailed  explanation  of  the  source  of  all  other  parameters  in  the  human  thrombopoiesis 
model. 


4.3.1  Maximum  mitotic  repopulation  rate  (a) 

The  maximum  mitotic  repopulation  rate  was  set  to  2.0  d^1  which  results  in  a  minimum 
division  time  of 


_  Z ri(2)  24 h 
~  2d-1  '  ~ld 


8.3  h 


(37) 


In  adult  humans,  cells  can  divide  up  to  about  twice  a  day,  and  the  DNA  synthesis  phase 
only  takes  approximately  8  h  (Andreeff  et  al.  2000).  Thus,  a  max  repopulation  rate  of  2.0 
d-1  is  biologically  reasonable. 


4.3.2  Platelet  transit  rate  (i/j) 

Several  studies  have  reported  a  platelet  lifespan  in  humans  of  approximately  9  d  (Schmitt 
et  ah  2001;  Mezzano  et  al.  1982;  Kaushansky  2005;  Finch  et  al.  1977).  Thus,  the  value  of  i/j 
was  set  to  1/9  or  0.11  d-1. 

4.3.3  MK  subcompartment  number  (n) 

The  number  of  megakaryocyte  subcompartments  was  set  to  10  (■ n  =  10).  Following  low 
dose  radiation  exposure  in  humans,  platelet  counts  remain  steady  for  a  period  of  time  before 
rapidly  declining.  These  results  can  be  interpreted  as  meaning  MKs  have  a  fairly  fixed 
maturation  time  with  little  variance,  suggesting  a  large  value  for  77.  However,  in  order  to 
decrease  computational  time  for  the  optimizations,  10  sub  compartments  was  chosen  as  an 
upper  threshold.  Increasing  n  beyond  10  only  improved  fits  slightly  but  increased  the  amount 
of  time  it  took  to  run  simulations  significantly. 


4.3.4  Platelet  subcompartment  number  (777) 

Fitting  human  platelet  data  to  gamma  functions  has  led  to  m  values  between  7  and  11  (Dowl¬ 
ing  et  al.  2010).  Here  777  was  set  to  the  midpoint  of  this  range,  9.  A  second  biomathematical 
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model  of  thrombopoiesis  used  7  platelet  subcompartments  (Scholz  et  al.  2010),  reinforcing 
the  notion  that  9  is  a  feasible  number. 

4.3.5  Ratio  of  receptor  concentrations  (03) 

The  relative  feedback  effect  of  x2  and  x3  on  the  x\  repopulation  rate  is  determined  by  the 
TPO  receptor  concentration  on  the  surface  of  MKs  and  platelets.  The  ratio  of  platelet  TPO 
receptors  to  MK  TPO  receptors  gives  an  approximate  value  of  the  ratio  0:>>  j ()2 .  In  humans, 
the  calculated  number  of  TPO  receptors  on  MKs  in  the  peripheral  and  umbilical  cord  blood 
is  1933  ±  772  (n=3)  and  184  ±  48  (n=4)  sites  per  cell,  respectively  (Kuwaki  et  al.  1998).  The 
average  value  of  the  peripheral  and  umbilical  cord  blood  MK  receptor  count  is  934  receptors 
per  MK.  The  number  of  TPO  receptors  in  healthy  humans,  using  data  from  5  experiments, 
was  30  ±  4  (Broudy  et  al.  1997).  Thus,  the  ratio  of  03  to  02  is  30/934  and  03  =  0.03  •  02. 

4.3.6  Maximum  and  minimum  rates  of  MK  maturation  ( Smin ,  Smax) 

The  maximum  and  minimum  rates  of  MK  maturation  were  determined  by  allowing  the 
immature  MK  maturation  time  to  decrease  or  increase  by  up  to  50%.  Thus,  if  the  immature 
MK  maturation  time  is  l/2<5o,  the  maximum  maturation  time  is  3/4<5o  and  the  minimum 
maturation  time  is  l/45o-  This  leads  to  minimum  and  maximum  total  MK  maturation  times 
of: 


3  1  _  5 

450  2<50  4<50 

1  1  _  3 

4<50  2<50  4  <50 


Thus, 


45o 

X 

4<5o 

IT 


(38) 

(39) 


(40) 

(41) 

(42) 


When  92%  of  platelets  were  removed  from  the  circulation  of  rats,  the  MK  maturation  time 
decreased  to  62-85%  of  control  values  (Ebbe  et  al.  1968).  This  provides  some  validation  for 
setting  Tmin  equal  to  75%  of  normal. 


4.3.7  Number  of  platelets  produced  per  MK  (a) 

The  number  of  platelets  produced  per  MK  (3000  platelets/MK)  was  assumed  to  be  the  same 
as  in  Smirnova’s  murine  and  human  models  (Smirnova  2011).  This  value  is  in  agreement 
with  other  studies  in  humans  (Deutsch  and  Tomer  2013). 
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4.3.8  Radiosensitivity  of  mitotic  progenitor  cells  (D®,  rii) 

In  humans,  the  Dq  and  n  values  for  MK  colony-forming  units  (CFU-MKs)  when  stimulated 
with  only  TPO  were  0.79  Gy  and  1.0,  respectively  (Kashiwakura  et  al.  2000).  These  values 
correspond  with  the  variables  for  mitotic  progenitor  radiosensitivity.  Thus,  -D^O.79  Gy  and 
ni=1.0  for  humans. 

4.3.9  Number  of  targets  per  MK  (n2) 

The  number  of  hits  required  to  kill  a  MK  was  set  to  4  (n2  =  4) .  This  variable  was  originally 
included  in  the  optimization  procedure,  but  there  was  a  large  collinearity  between  ?72  and 
D®  making  it  not  possible  to  optimize  these  two  parameters  together.  Thus,  rather  than 
optimize  ?72,  the  value  was  systematically  increased  until  a  decline  was  seen  in  the  early 
plateau  phase  at  a  dose  level  supported  by  experimental  data. 

4.3.10  Effect  of  damaged  states  on  mediator  levels  (fi) 

The  death  rate  of  damaged  cells  was  set  to  1.0  d-1.  In  Smirnova’s  original  model,  the 
moderately  and  heavily  damaged  compartments  had  death  rates  of  0.5  d-1  and  6.0  d^1  to 
simulate  mitotic  and  interphase  death,  respectively  (Smirnova  2011).  Here,  the  damaged 
compartment  contains  both  cells  that  die  during  mitosis  and  interphase.  Thus,  the  value  of 
1.0  d-1  provides  an  approximate  average  lifespan  of  1.0  d. 

4.4  Optimization  Data 

For  the  radiation  model  optimization,  data  was  collected  from  multiple  sources.  Table  4.2 
provides  a  summary  of  the  radiation  data  used  in  optimizing  the  human  thrombopoiesis 
model.  Information  on  the  data  reference,  dose,  and  baseline  values  is  given  for  all  case 
studies.  Data  from  17  case  studies  were  used  in  the  optimization,  11  of  which  contained 
data  points  up  to  and  beyond  30  days  post-radiation  exposure.  The  subjects  received  doses 
ranging  from  0.12  to  12.5  Gy. 

4.5  Results:  Optimization  and  MCMC  Analysis 

The  optimized  parameter  values  and  corresponding  95%  confidence  intervals  are  given  in 
Table  4.3.  Also  provided  are  the  initialization  value,  the  lower  bound,  and  the  upper  bound 
placed  on  each  parameter  prior  to  running  the  optimization.  Table  4.4  gives  a  complete 
list  of  parameter  values,  descriptions,  and  bases.  Figure  4.3  shows  the  optimized  model 
compared  to  select  data  used  for  the  optimization.  Each  dataset  comes  from  the  specified 
case  study/patient.  The  remainder  of  the  optimization  data  is  compared  to  the  model  output 
in  Appendix  B.l.  For  details  on  the  MCMC  results  and  identifiability  analysis  see  Appendix 
C. 
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Table  4.2:  Optimization  data  summary  for  thrombopoiesis  model 


Incident /Subject 

Reference 

Dose  (Gy) 

Baseline  (10 3/iL  x) 

Los  Alamos  1945-2 

Hempelmann  et  ah  1952 

0.12 

819.26“ 

UT  CARL  1971 

Andrews  et  ah  1961* 

2.6 

350.63“ 

Argonne  1952-1 

Bond  et  ah  1965 

1.5 

111.23“ 

Y-12  1958-B 

Andrews  et  ah  1961* 

2.97 

161.25“ 

Y-12  1958-C 

Andrews  et  ah  1961* 

3.73 

181.25“ 

Y-12  1958-A 

Bond  et  ah  1965 

4.02 

289.56“ 

Vinca  1958-B 

Jarnmet  et  ah  1959 

3.51 

210.31“ 

Vinca  1958-H 

Jarnmet  et  ah  1959 

4.37 

233.73“ 

Vinca  1958-M 

Jarnmet  et  ah  1959 

5.74 

219.88“ 

Vinca  1958-G 

Jarnmet  et  ah  1959 

5.35 

250. 00b 

Vinca  1958-D 

Jarnmet  et  ah  1959 

5.40 

250. 006 

Vinca  1958-V 

Jarnmet  et  ah  1959 

6.11 

221.29“ 

Yarmonenko  1988-Z 

Yarmonenko  1988 

9.8 

140.21“ 

Nesvizh  1991 

Baranov  et  ah  1994 

12.5 

106.75“ 

Cancer  patients  1958- 
mean  of  18 

Miller  et  ah  1958 

1.0 

197.25c 

Cancer  patients  1958- 
mean  of  12 

Miller  et  ah  1958 

1.5 

234. 25c 

Cancer  patients  1958- 
mean  of  30 

Miller  et  ah  1958 

2.0 

155. 75c 

“Mean  of  up  to  four  data  points  collected  within  3  days  of  radiation  exposure 
6Mean  platelet  concentration  observed  in  humans  (Valentin  2002) 

“Mean  of  pre-radiation  data  points 
*Data  provided  by  Ron  Goans 

Additional  details  on  the  human  case  studies  provided  in  Appendix  A. 

4.6  Validation 

To  validate  the  human  thrombopoiesis  model,  the  model  outputs  were  compared  to  data 
on  platelet  count  following  radiation  in  humans  that  were  not  included  in  the  optimization 
(Hempelmann  et  ah  1952;  Bond  et  al.  1965;  Hirama  et  ah  2003;  Stavem  et  ah  1985;  Mettler 
et  ah  2001;  Mettler  2001).  Select  validation  results  are  shown  in  Figure  4.4.  The  remainder 
of  the  validation  results  are  presented  in  Appendix  B.2. 

Validation  was  also  performed  by  comparing  model  output  to  kinetic  platelet  data  from 
patients  following  thrombopheresis  (Lasky  et  ah  1981;  Weisbach  et  ah  1999;  Sullivan  et  ah 
1977).  In  thrombopheresis  a  certain  percentage  of  platelets  are  removed  from  the  blood 
stream.  We  can  initialize  the  value  of  x 3  in  our  model  using  this  percentage,  where  we 
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Table  4.3:  Thrombopoiesis  optimized  parameters  and  confidence  intervals 


Parameter  Optimized  Value 

2.5%“ 

97. 5%6 

Initial  Value  (Lower, Upper)c 

7  (d-1) 

0.2556 

0.2254 

0.2819 

0.2647 

(0.01,  1.99) 

6  (d-1) 

0.0589 

0.0503 

0.0683 

0.0586 

(0.05,  0.2) 

e2 

0.0063 

0.0041 

0.0109 

0.0049 

(0.0039,  20) 

X 

3.51 

3.4772 

3.5293 

3.4913 

(0.3,  4) 

D°2  (Gy) 

2.768 

2.5914 

2.9329 

2.8794 

(1,  50) 

“Lower  end  of  parameter  95%  confidence  interval 
bUpper  end  of  parameter  95%  confidence  interval 

“Initialization  value,  lower  (min),  and  upper  (max)  bound  given  to  parameter  during  optimization 


assume  (1)  an  equal  fraction  of  platelets  is  removed  from  each  age  dependent  platelet  sub¬ 
compartment  and  (2)  the  removal  of  platelets  is  instantaneous.  In  one  study,  thrombopheresis 
was  performed  twice.  To  simulate  this,  we  initialized  the  value  of  according  to  the  first 
thrombopheresis  and  manually  altered  Xs  at  the  time  of  the  second  thrombopheresis  without 
changing  the  values  of  the  other  state  variables.  This  validation  helps  to  confirm  that  the 
kinetic  rates  of  cell  maturation  and  feedback  are  correct  in  the  model.  These  results  are 
shown  in  Figure  4.5. 

4.6.1  Relative  compartment  sizes 

To  further  validate  the  thrombopoiesis  model,  the  relative  compartment  sizes  at  equilibrium 
are  calculated  using  model  parameters  and  compared  to  experimental  data.  An  equilibrium 
point  of  the  model  is  a  point  in  the  steady  state  at  which  the  vector  field  vanishes.  Therefore, 
Eq.  25-30  imply  that  at  an  equilibrium  point 

0  =  -  nSx^A 

0  =  nSx^j  ^i  —  nSx^j  i  ( i  =  2,  3, ...,  n/2) 

0  =  nSx^n/2  -  tiSqX 

0  =  n5 0Z2M,i-i  -  n80x% M,i  (*  =  2, 3, ...,  n/2) 

0  =  anS0x 2M,n/2  -  mi>x 3,1 

0  =  m^x^i-i  —  mil)Xzti  ( i  =  2, 3, ...,  m) 

Moreover,  at  equilibrium  S  =  S0,  so  we  must  have 

x™d  nS0 
*21,1  7 

r£U( i 

4^  =  1  (i  =  1,  2, ...,  n/2;  j  =  1,  2, ...,  n/2) 

X2I,j 
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Table  4.4:  Biological  parameter  values  for  thrombopoiesis 


Parameter 

Description 

Value 

Section  (References) 

a 

Maximum  rate  of  mitotic  progenitor 
cell  repopulation 

2.0  d-1 

4.3.1 

7 

Rate  of  mitotic  progenitor  cell  matu¬ 
ration 

0.26  d”1 

Optimized 

<5o 

Rate  of  MK  maturation 

0.06  d"1 

Optimized 

Rate  of  platelet  decay 

0.11  d-1 

4.3.2  (Mezzano  et  al.  1982; 
Finch  et  al.  1977) 

n 

Number  of  MK  subcompartments 

10 

4.3.3 

m 

Number  of  platelet  subcompartments 

9 

4.3.4  (Dowling  et  al.  2010; 
Scholz  et  al.  2010) 

e2 

Decay  rate  of  mediator  due  to  x2  cells 
relative  to  decay  rate  due  to  X\  cells 

0.006 

Optimized 

6s 

Decay  rate  of  mediator  due  to  x%  cells 
relative  by  decay  rate  due  to  x±  cells 

0.00018 

4.3.5  (Kuwaki  et  al.  1998; 
Broudy  et  al.  1997) 

& min 

Minimum  rate  of  MK  maturation 

0.048 

4.3.6 

$max 

Maximum  rate  of  MK  maturation 

0.08 

4.3.6 

A 

Strength  of  X3  feedback  on  immature 
MK  maturation 

3.5 

Optimized 

a 

Number  of  platelets  per  MK 

3000 

4.3.7  (Smirnova  2011; 
Deutsch  and  Tomer  2013) 

D ? 

Determines  fraction  of  X\  targets  that 
are  hit  by  radiation  at  dose  D 

0.79  Gy 

4.3.8  (Kashiwakura  et  al. 
2000) 

ni 

Number  of  targets  per  x\  cell 

1 

4.3.8  (Kashiwakura  et  al. 
2000) 

n° 

u2 

Determines  fraction  of  x2  targets  that 
are  hit  by  radiation  at  dose  D 

2.77  Gy 

Optimized 

n2 

Number  of  targets  per  x2  cell 

4 

4.3.9 

h 

Rate  of  moderately  damaged  cell 
death 

1.0  d-1 

4.3.10 

23 


c/5 

c/5 

05 

c= 

o 

C/5 

c= 

0) 

E 

33 

c= 

o 

-I — 1 

cc! 

i— 
H — ' 

c= 

CD 

O 

c= 

o 

o 

-I — ' 

Q 

CD 

-t — > 
0- 


(a)  Dose:  1.5  Gy  (1-2  Gy;  Argonnne) 


(b)  Dose:  2.97  Gy  ±  25% 


(c)  Dose:  4.02  Gy  ±  25% 


(d)  Dose:  9.8  Gy 


Figure  4.3:  Thrombopoiesis  model  compared  to  select  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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(a)  Dose:  1.68  Gy 


(b)  Dose:  2.59  Gy  ±  25% 


(c)  Dose:  3.6  Gy  (d)  Dose:  12.0  Gy 

Figure  4.4:  Thrombopoiesis  model  compared  to  select  validation  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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Platelet  Concentration  (dimensionless)  _  Platelet  Concentration  (dimensionless) 


a)  Lasky  et  at.  1981;  multiple  subjects  (b)  Weisbach  et  al.  1999;  22  subjects 


(c)  Sullivan  et  al.  1977;  individual  data  (d)  Sullivan  et  al.  1977;  individual  data 

Figure  4.5:  Thrombopoiesis  model  overlaid  on  thrombopheresis  data 
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the  relative  compartment  sizes  at  equilibrium  are 
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In  humans,  the  total  number  of  bone  marrow  cells  is  approximately  2  •  1010  per  kg  of 
body  weight  (Valentin  2002),  and  the  fraction  of  those  estimated  to  be  MKs  is  0.04%  (Levine 
1980).  This  gives  a  total  MK  count  of  8  •  106  per  kg  body  weight.  The  number  of  platelets 
circulating  is  approximately  250  •  103  per  /i L  of  blood  (Valentin  2002).  Humans  have  71  mL 
of  blood  per  kg  of  body  weight  (Starr  and  Taggart  1989),  leading  to  a  total  platelet  count 
of  1.775  •  1010  cells  per  kg  of  body  weight.  This  leads  to  an  X2  to  x3  ratio  of  0.0005,  which 
is  comparable  with  the  model  prediction  of  0.0006. 
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5  Granulopoiesis  Model 

5.1  Conceptual  Model 

Granulopoiesis  is  the  process  by  which  mature  granulocytes  are  generated  from  pluripo- 
tent  HSCs.  Granulocytes  are  a  type  of  leukocyte  and  can  be  subdivided  into  neutrophils, 
eosinophils,  and  basophils.  Neutrophils  account  for  the  majority  of  granulocytes  and  are 
involved  in  phagocytosis,  the  release  of  soluble  anti-microbials,  and  generating  neutrophil 
extracellular  traps  (NETs). 

The  granulopoiesis  model  has  four  compartments:  mitotic  precursors  in  the  bone  marrow 
(Xj ),  post-mitotic  precursors  in  the  bone  marrow  (X2),  granulocytes  in  circulation  (X3),  and 
granulocytes  in  the  tissues  (A"4).  Figure  5.1  illustrates  the  structure  of  the  granulopoiesis 
model.  In  this  diagram,  the  damaged  compartments  represent  both  the  weakly  damaged 
and  the  damaged  cells  in  the  system. 


Cell  transitions 


-  ■  Instant  radiation  damage 


■>  Feedback  effects 


Figure  5.1:  Granulopoiesis  model  diagram 


5.1.1  Radiation  damage 

All  cells  in  the  granulopoietic  lineage  are  radiosensitive.  The  mitotic  cells  in  the  bone  marrow 
are  the  most  radiosensitive  and  become  either  damaged  or  weakly  damaged.  The  cells  in 
the  other  three  compartments  are  less  radiosensitive  than  mitotic  cells  and  enter  a  single 
damaged  state  following  radiation  exposure. 


29 


The  weakly  damaged  cell  compartment  is  used  to  model  the  late  transient  increase  in 
granulocyte  counts  observed  following  radiation  exposure  (Bond  et  al.  1965).  Weakly  dam¬ 
aged  cells  reproduce  like  normal  cells  for  a  period  of  time  defined  by  the  following  equation: 

A tae  =  r  -vD  (43) 


where  D  is  the  radiation  dose,  and  r  and  v  are  parameters.  Once  the  elapsed  time  after 
radiation  exposure  is  equal  to  A tae,  all  weakly  damaged  cells  and  their  progeny  start  decaying 
at  rate  r/. 

5.1.2  Feedback  mechanisms 

The  model  has  two  feedback  mechanisms.  The  first  involves  the  generic  mediator  which  stim¬ 
ulates  the  repopulation  rate  of  X{ .  The  mediator’s  concentration  decreases  as  the  population 
of  granulopoietic  cells  increases.  Because  the  action  of  the  mediator  is  fast,  its  concentration 
can  be  written  as  a  function  of  cell  concentration,  and  thus  the  rate  of  repopulation  (B)  can 
also  be  written  as  a  function  of  cell  concentration  (see  Section  3.2  for  derivation): 


~  1  +  P(xf  +  xf  +  x™d  +  Y^t=2ei(xid  +  xt  +  xid)) 

The  second  feedback  mechanism  affects  the  rate  of  granulocyte  transition  from  the  bone 
marrow  into  the  blood.  As  the  concentration  of  granulocytes  in  the  blood  increases,  the  rate 
of  transition  from  the  bone  marrow  (F)  decreases  resulting  in  a  negative  feedback  loop.  The 
equation  for  F  is 


p  =  jl  +  M(x¥  +  xi  +  x%d)2 

1  +  L(xf  +  x$  +  x%d)2  {  J 

In  the  model,  rather  than  defining  M  and  L  to  describe  the  granulocyte  rate  of  entry  from 
the  bone  marrow  to  the  blood  (F),  m  and  l  are  used,  m  is  equivalent  to  M{x^d)2  and  l  is 
equivalent  to  L(x 3d)2  where  X3  represents  the  concentration  of  X3  at  equilibrium. 

The  sources  for  parameters  defining  the  rates  between  compartments  (7,  S,  m,  l,  k),  the 
rates  of  cell  decay  (^,  /i,  7),  radiosensitivity  parameters  (D®,  D °00,  rii,  D n2,  D®,  n3,  D®, 
n4),  and  other  parameters  ( a ,  d2,  03,  64,  r,  v )  are  given  in  Section  5.3. 

5.2  Mathematical  Model 

The  mathematical  model  of  granulopoiesis  is  described  by  a  set  of  initial  conditions  and  dif¬ 
ferential  equations  that  describe  the  change  in  cell  concentrations  over  time.  The  description 
of  each  cell  type  is  given  in  Table  5.1. 

At  equilibrium  cell  concentrations  in  the  Xi,  X2,  X3,  and  X4  compartment  are  equal  to 
X}d,  x%d,  x%d,  and  respectively.  Following  acute  exposure,  cells  immediately  transition 
into  the  damaged  or  weakly  damaged  compartments  depending  on  their  radiosensitivity. 
$  determines  the  ratio  of  the  number  of  cells  that  transition  into  the  damaged  or  weakly 
damaged  compartment.  The  following  equations  describe  how  cell  concentrations  change 
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Table  5.1:  Cell  concentrations  in  each  compartment  of  the  granulopoiesis  model 


Variable  name  Description 

Xid  Undamaged  cells  in  the  mitotic  compartment  X\ 

x1^1  Undamaged  cells  in  the  post  mitotic  compartment  X2 

x^d  Undamaged  granulocytes  in  the  blood  A"3 

xvAd  Undamaged  granulocytes  in  the  tissues  A"4 

xf  Damaged  cells  in  the  mitotic  compartment  Ad 

xd  Damaged  cells  in  the  post  mitotic  compartment  X2 

xd  Damaged  granulocytes  in  the  blood  A3 

x4  Damaged  granulocytes  in  the  tissue  X4 

xfd  Weakly  damaged  cells  in  the  mitotic  compartment  A4 

X2d  Weakly  damaged  cells  in  the  post  mitotic  compartment  X2 

x™d  Weakly  damaged  granulocytes  in  the  blood  A"3 

x™d  Weakly  damaged  granulocytes  in  the  tissue  X4 

x^d  Concentration  of  cells  in  the  A,;  compartment  at  equilibrium  (i=l,2,3,4) 


following  an  acute  exposure  of  dose  D  and  determine  initial  conditions  for  the  temporal 
model. 


where 


f  =  xf  (l  -  (1  -  e~D/D>)ni )  (i  =  1,  2,  3, 4) 
d  (l  -e-D/Dty1 


xwd  _  x ud 


1  ~  1  1+7? 
x?d  =  0  (i  =  2,3,4) 
1 


W  _  -ud 

Jb  ^  JU  ^ 


1  +t9 


(l-e-wof)' 


(*  =  2,3,4) 


0  = 


(1  -  e_D''Df)“1  -  (1  -  e_D/D“0)”‘ 
(1  _  e-D/D?“)m 


(46) 

(47) 

(48) 

(49) 

(50) 


Once  radiation  exposure  occurs,  a  series  of  differential  equations  are  used  to  describe 
the  compartmental  transitions,  repopulation  and  decay  rates,  and  the  feedback  mechanisms 
shown  in  Figure  5.1.  The  following  equations  describe  the  model  where  Eq.  46-50  are  used 
to  calculate  initial  conditions,  iy  represents  dxi/dt. 


x?  =  Bxf  -  7xf 


(51) 
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(52) 

xf  =  Fxf  -  nxf 

(53) 

-rud  —  KTud  —  'ihrrud 

Jb  4  —  rXjJbr^  (/y  0C4 

(54) 

±f  =  j 

\  Bxf  -  'yxf 
—T)xf 

to  <  t  <  to  +  Afae 

t>tQ  +  Atae 

(55) 

if  =j 

\  ^xf  —  Fxf 

1  —T]xf 

t0  <  t  <  t0  +  Atae 
t  >  t0  +  Atae 

(56) 

if  =| 

\  Fxf  —  nxf 

1  —T]xf 

tQ  <  t  <  t0  +  Atae 
t>tQ  +  Atae 

(57) 

if  =| 

\  K,xf  —  pixf 

^  —rjxf 

t0<t<t0  +  Atae 
t>tQ  +  Atae 

(58) 

xd  =  - 

-/ixf  (i  —  1,  2, 

3,4) 

(59) 

where 


B 

F 

M 

L 


a 

1  +  /3(xf  +  xf  +  xf  +  J2i=:  2  6i(xid  +  Xd  +  x™d)) 

1  +  M  (xf  +  xf  +  x™  d)2 

1  +  L(x3d  +  xf  +  X3 ' d)2 
m 


(xf)2 

l 


(xf)2 


(60) 

(61) 

(62) 

(63) 


5.3  Radiation  Model  Parameter  Sources 

Granulopoiesis  parameters  were  determined  using  experimental  data  or  through  optimiza¬ 
tion.  The  optimized  parameters  include  ip,  l,  D®,  D((00,  D®  and  77.  The  following  gives  a 
detailed  explanation  of  the  source  of  all  other  parameters  in  the  human  granulopoiesis  model. 


5.3.1  Maximum  mitotic  repopulation  rate  (a) 

The  maximum  mitotic  repopulation  rate  was  set  to  2.0  d^1  which  results  in  a  minimum 
division  time  of 


_  ln( 2)  24h 
“  2d-1  '  ~ld 


8.3  h 


(64) 


In  adult  humans,  cells  can  divide  up  to  about  twice  a  day,  and  the  DNA  synthesis  phase 
only  takes  approximately  8  h  (Andreeff  et  al.  2000).  Thus,  a  max  repopulation  rate  of  2.0 
d_1  is  biologically  reasonable. 
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5.3.2  Compartment  transition  rates  (6,  m,  l,  k,  7) 

The  mitotic  and  post-mitotic  transition  times  of  granulocytes  have  been  experimentally 
determined  (Price  et  al.  1996).  The  mitotic  transition  time  is  approximately  6.0  d,  leading 
to  7  =  1/6  =  0.17  d-1.  The  post-mitotic  transition  time  is  5.5  d  leading  to  heg=l/6.4=0.16 
d-1.  Following  G-CSF  treatment,  the  minimum  transit  time  is  2.9  d,  leading  to  a  maximum 
transition  rate,  or  5,  of  0.34  d-1. 

Because 


^  =  5^  (65) 

we  can  set  m  as  a  function  of  Seq,  <5,  and  l. 

m  =  %(  1  +  0-1  (66) 

0 

5eq  and  5  have  already  been  determined,  l  was  included  in  the  optimization  procedure,  and 
m  was  set  using  the  above  equation.  The  granulocyte  half-life  in  the  blood  of  humans  is  10.4 
h  leading  to  a  transit  rate  n  of  1.6  d-1  (Price  et  al.  1996). 

5.3.3  Receptor  concentration  ratios  ( 02 ,  #3,  d4) 

The  ratio  of  receptors  on  x2  and  X3  cells  compared  with  X\  cells  was  determined  based 
on  experimental  data  on  the  number  of  G-CSF  receptors  per  granulocyte  progenitor  and 
peripheral  granulocyte  (Shinjo  et  al.  1995).  Granulocyte  progenitors  are  known  to  express 
both  the  CD33  and  CD34  antigens.  Using  these  markers  for  granulocyte  progenitor  selection, 
the  relative  number  of  G-CSF  receptors  on  early  progenitors  in  humans  was  measured  to 
be  0.0446  ±  0.0068  and  on  later  progenitors  was  0.0988  ±  0.0250.  The  relative  number  of 
G-CSF  receptors  on  peripheral  granulocytes  in  humans  was  8353  ±  623.  This  gives  ratios 
of  998/446=2.24  and  8353/988=8.45  for  02  and  #3,  respectively.  Because  A"3  and  X4  both 
represent  granulocytes,  04  was  assumed  to  equal  #3. 

5.3.4  Granulocyte  radiosensitivity  (D 3,  D4) 

Granulocyte  radiosensitivity  was  assumed  to  be  equal  to  the  post-mitotic  cell  radiosensitivity 
(H3  =  D4  =  D °).  In  Smirnova’s  original  model,  all  three  cell  types  also  have  the  same 
radiosensitivity  (Smirnova  2011). 

5.3.5  Time  of  abortive  rise  Atae  (r,  v ) 

The  time  of  abortive  elevation  A tae  is  determined  using  the  following  equation: 

Atae  =  t  vD  (67) 


where  D  is  the  radiation  dose,  r  and  v  were  determined  through  examination  of  the  time  of 
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abortive  rise  in  the  granulopoiesis  optimization  data  (see  Section  5.4).  In  each  case  study, 
an  approximate  time  of  abortive  rise  was  graphed  against  the  acute  radiation  dose.  Fitting 
Equation  67  to  the  data,  resulted  in  values  of  v  =  1.25  d  Gy-1  and  r  =  22.70  d. 

5.3.6  Targets  per  cell  (n\,  n-2,  n^,  n^) 

The  number  of  targets  per  cell  in  the  granulopoiesis  model  was  assumed  to  be  1  for  all 
cells  in  the  lineage  (nj=l  for  %  =  1,2,  3, 4)  .  Using  this  assumption,  we  are  simplifying  our 
radiation  model  to  the  single-target,  single-hit  model  used  by  Smirnova  (Smirnova  2011). 
Using  this  model  of  radiation  damage,  we  are  able  to  accurately  match  cell  survivability 
following  radiation  exposure. 

5.3.7  Effect  of  damaged  states  on  mediator  levels  (fi) 

The  death  rate  of  damaged  cells  was  set  to  1.0  d_1.  In  Smirnova’s  original  model,  the 
moderately  and  heavily  damaged  compartments  had  death  rates  of  0.5  d-1  and  6.0  d^1  to 
simulate  mitotic  and  interphase  death,  respectively  (Smirnova  2011).  Here,  the  damaged 
compartment  contains  both  cells  that  die  during  mitosis  and  interphase.  Thus,  the  value  of 
1.0  d-1  provides  an  approximate  average  lifespan  of  1.0  d. 

5.4  Optimization  Data 

For  the  radiation  model  optimization,  the  data  was  collected  from  multiple  sources.  Table 
5.2  provides  a  summary  of  the  radiation  data  used  in  optimizing  the  human  granulopoiesis 
model.  Information  on  the  data  reference,  dose,  and  baseline  values  is  given  for  all  case 
studies.  Data  from  16  case  studies  were  used  in  the  optimization,  11  of  which  contained 
data  points  up  to  and  beyond  30  days  post-radiation  exposure.  The  subjects  received  doses 
ranging  from  0.12  to  9.8  Gy. 

5.5  Results:  Optimization  and  MCMC  Analysis 

The  optimized  parameter  values  and  corresponding  95%  confidence  intervals  are  given  in 
Table  5.3.  Also  provided  are  the  initialization  value,  the  lower  bound,  and  the  upper  bound 
placed  on  each  parameter  prior  to  running  the  optimization.  Table  5.4  gives  a  complete 
list  of  parameter  values,  descriptions,  and  bases.  Figure  5.2  shows  the  optimized  model 
compared  to  select  data  used  for  the  optimization.  Each  dataset  comes  from  the  specified 
case  study /patient.  The  remainder  of  the  optimization  data  is  compared  to  the  model 
output  either  in  the  Discussion  (Section  7;  Figure  7.3b)  or  in  Appendix  B.l.  For  details  on 
the  MCMC  results  and  identihability  analysis  see  Appendix  C. 

5.6  Validation 

To  validate  the  human  granulopoiesis  model,  the  model  outputs  were  compared  to  data  on 
granulocyte  count  following  radiation  in  humans  that  were  not  included  in  the  optimization 
(Hempelmann  et  al.  1952;  Bond  et  al.  1965;  Hirama  et  al.  2003;  Stavern  et  ah  1985;  Mettler 
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Table  5.2:  Optimization  data  summary  for  granulopoiesis  model 


Incident /Subject 

Reference 

Dose  (Gy) 

Baseline  (103/iL  3) 

Los  Alamos  1945-2 

Hempelmann  et  al.  1952 

0.12 

4.102a 

UT  CARL  1971 

Andrews  et  al.  1961* 

2.6 

4.102a 

Argonne  1952-1 

Bond  et  al.  1965 

1.5 

4.102a 

Y-12  1958-B 

Andrews  et  al.  1961* 

2.97 

4.102a 

Y-12  1958-C 

Andrews  et  al.  1961* 

3.73 

4.102a 

Y-12  1958-A 

Bond  et  al.  1965 

4.02 

4.102a 

Vinca  1958-B 

Jammet  et  al.  1959 

3.51 

4.102a 

Vinca  1958-H 

Jammet  et  al.  1959 

4.37 

4.102a 

Vinca  1958-M 

Jammet  et  al.  1959 

5.74 

4.102a 

Vinca  1958-G 

Jammet  et  al.  1959 

5.35 

4.102a 

Vinca  1958-D 

Jammet  et  al.  1959 

5.40 

4.102a 

Vinca  1958-V 

Jammet  et  al.  1959 

3.11 

4.102a 

Yarmonenko  1988-Z 

Yarmonenko  1988 

9.8 

4.102a 

Cancer  patients  1958- 
mean  of  18 

Miller  et  al.  1958 

1.0 

6.22b 

Cancer  patients  1958- 
mean  of  12 

Miller  et  al.  1958 

1.5 

6.99b 

Cancer  patients  1958- 
mean  of  30 

Miller  et  al.  1958 

2.0 

8.13b 

“Mean  granulocyte  concentration  observed  in  humans  (Valentin  2002) 
bMean  of  pre-radiation  data  points 
*Data  provided  by  Ron  Goans 

Additional  details  on  the  human  case  studies  provided  in  Appendix  A 

et  al.  2001;  Mettler  2001).  Select  validation  results  are  shown  in  Figure  5.3.  The  remainder 
of  the  validation  results  are  given  in  Appendix  B.2. 

5.6.1  Relative  compartment  sizes 

To  further  validate  the  granulopoiesis  model,  the  relative  compartment  sizes  at  equilibrium 
are  calculated  using  model  parameters  and  compared  to  experimental  data.  An  equilibrium 
point  of  the  model  is  a  point  in  the  state  space  at  which  the  vector  field  vanishes.  Therefore, 
Eq.  52  -  54  imply  that  at  an  equilibrium  point 

0  =  -  6- 

0  =  .) '  +  £?  -  Kxf 

l+l  2  3 
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Table  5.3:  Granulopoiesis  optimized  parameters  and  confidence  intervals 


Parameter 

Optimized  Value 

2.5%“ 

97.5%fc 

Initial  Value  (Lower, Upper)c 

^  (d 

-1) 

1.0542 

0.4409 

1.6179 

11.8422  (0.01,  20) 

l 

1.6906 

1.1894 

2.3066 

1.1539  (1.126,  12) 

D000 

(Gy) 

2.4591 

1.9941 

2.7815 

2.3881  (1,  10) 

D°2  (Gy) 

45.76 

43.3506 

47.2832 

10  (2,  50) 

V  (d~ 

-1) 

0.08161 

0.0524 

0.1645 

0.0865  (0.05,  5) 

D°,  (Gy) 

0.4605 

0.3093 

0.6766 

0.4714  (0.3,  0.9) 

“Lower  end  of  parameter  95%  confidence  interval 
&Upper  end  of  parameter  95%  confidence  interval 

“Initialization  value,  lower  (min),  and  upper  (max)  bound  given  to  parameter  during  optimization 


0  =  KX^d  —  tpX^1' 


Therefore,  at  equilibrium  we  must  have 


~ud 

1 

~ud 

x2 


X  1 +m 
U  l+l 

7 

0.34{±^ 

1+1’69  =  0.94 


0.17 


pz.ud 

X2 


K 


X 


ud 


:  1+m 


1  -\-l 

1-60 

— ^  - ii  n  T7  —  9.97 


~ud 

J/4 


x 


ud 


— > 


0.34 

K 

1.60 

L05 


1+0.27 

1+1.69 


=  1.52 


In  human  granulopoiesis  the  total  number  of  mitotic  cell  progenitors  is  2.11  •  109  per  kg  of 
body  weight,  the  total  number  of  post-mitotic  progenitors  is  5.59  •  109  per  kg  of  body  weight, 
and  the  total  number  of  granulocytes  in  the  blood  (both  marginated  and  circulating)  is 
0.61  •  109  per  kg  of  body  weight  (Dancey  et  al.  1976).  This  leads  to  the  following  ratios: 


Xi 

x2 

Xi 

x3 


0.377 

3.459 


(68) 

(69) 
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X2 

X3 


9.164 


(70) 


Based  on  this  experimental  data,  the  ratio  given  by  the  kinetic  parameters  of  the  gran¬ 
ulopoiesis  model  defining  x2  to  x3  is  approximately  correct.  The  ratio  of  x\  to  x2  predicted 
by  the  model  is  larger  than  the  ratio  predicted  by  the  experimental  data,  suggesting  the 
value  of  7  used  by  the  model  is  too  small.  However,  the  value  of  7  used  in  our  model  has 
support  from  the  literature  (Price  et  al.  1996).  The  mitotic  cell  compartment  in  our  model 
also  contains  early  HSC  progenitors,  which  may  not  be  included  in  the  granulocyte  mitotic 
cell  progenitor  count  given  by  Dancey  et  al.  1976.  If  this  were  the  case,  the  calculated  ratio 
of  xi  to  x2,  based  on  relative  cell  compartment  sizes,  would  be  larger  than  0.377  and  closer 
to  agreement  with  the  model  prediction  of  0.94. 
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Table  5.4:  Biological  parameter  values  for  granulopoiesis 


Parameter  Description 

Value 

Section  (References) 

a 

Maximum  rate  of  mitotic  progenitor 

2.0  d”1 

5.3.1  (Andreeff  et  al. 

cell  repopulation 

2000) 

7 

Rate  of  mitotic  progenitor  cell  matura¬ 
tion 

0.17  d”1 

5.3.2  (Price  et  al.  1996) 

5 

Determines:  Post-mitotic  maturation 

rate 

0.34  d"1 

5.3.2  (Price  et  al.  1996) 

771 

Determines:  Post-mitotic  maturation 

rate 

0.27 

5.3.2  (Price  et  al.  1996) 

l 

Determines:  Post-mitotic  maturation 

rate 

1.69 

Optimized 

hi 

Rate  of  transition  from  blood  to  tissue 

1.60  d"1 

5.3.2  (Price  et  al.  1996) 

7\) 

Rate  of  granulocyte  decay  from  tissue 

1.05  d”1 

Optimized 

V 

Rate  of  weakly  damaged  cell  death 
(when  t  >  A tae) 

0.08  d”1 

Optimized 

e2 

Decay  rate  of  mediator  due  to  x2  cells 
relative  to  decay  rate  due  to  £1  cells 

2.24 

5.3.3  (Shinjo  et  al.  1995) 

o3 

Decay  rate  of  mediator  due  to  cells 

relative  to  decay  rate  due  to  x\  cells 

8.45 

5.3.3  (Shinjo  et  al.  1995) 

<h 

Decay  rate  of  mediator  due  to  x 4  cells 
relative  to  decay  rate  due  to  Xi  cells 

8.45 

5.3.3  (Shinjo  et  al.  1995) 

D\ 

Determines  fraction  of  damaged  X\  cells 

0.46  Gy 

Optimized 

Dooo 

Determines  ratio  of  weakly  damaged  to 
damaged  x\  cells. 

2.46  Gy 

Optimized 

D° 

Determines  fraction  of  damaged  x2  cells 

45.76  Gy 

Optimized 

n° 

-^3 

Determines  fraction  of  damaged  x%  cells 

45.76  Gy 

5.3.4 

n° 

u  4 

Determines  fraction  of  damaged  x 4  cells 

45.76  Gy 

5.3.4 

T 

Determines  time  of  abortive  rise: 

A tae  =  T  ~VD 

22.70  d 

5.3.5 

V 

Determines  time  of  abortive  rise: 

A tae  =  T  ~  VD 

1.25  d  Gy"1 

5.3.5 

rii 

Number  of  targets  per  cell  x% 

(z  =  1,2,3,4) 

1 

5.3.6 

n 

Rate  of  damaged  cell  death 

1.0  d”1 

5.3.7 
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CO 


co 


(a)  Dose:  1.5  Gy  (1-2  Gy;  Argonne) 


(b)  Dose:  2.97  Gy  ±  25% 


(c)  Dose:  3.51  Gy  ±  15%  (d)  Dose:  9.8  Gy 


Figure  5.2:  Granulopoiesis  model  compared  to  select  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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05 


Time  following  radiation  (days) 


(a)  Dose:  1.1  Gy 

05 


Time  following  radiation  (days) 


(c)  Dose:  3.6  Gy 


05 


(b)  Dose:  2.59  Gy  ±  25% 


05 


(d)  Dose:  12  Gy 


Figure  5.3:  Granulopoiesis  model  compared  to  select  validation  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region.  If  a  generic  baseline  value  was  used,  output  initialized 
at  the  upper  and  lower  end  of  the  granulocyte  range  is  shown  by  dotted  lines. 
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6  Lymphopoiesis  Model 

6.1  Conceptual  Model 

Lymphopoiesis  is  the  process  by  which  lymphocytes  are  generated  from  the  bone  marrow. 
The  three  main  types  of  lymphocytes  are  T  cells,  B  cells,  and  natural  killer  (NK)  cells.  B 
and  NK  cells  mature  in  the  bone  marrow,  while  T  cell  progenitors  migrate  to  the  thymus 
for  maturation.  Once  mature,  cells  enter  the  circulation  and  peripheral  lymphoid  organs 
including  the  spleen  and  lymph  nodes.  The  human  lymphopoiesis  model  consists  of  three 
compartments:  mitotic  precursors  in  the  bone  marrow  (Ad),  post-mitotic  precursors  in  the 
bone  marrow  or  thymus  (Ad),  and  lymphocytes  in  circulation  (Ad).  Figure  6.1  illustrates 
the  structure  of  the  model. 


- Cell  transitions  -  ■  ■»  instant  radiation  damage  . ■>  Feedback  effects 

Figure  6.1:  Lymphopoiesis  model  diagram 


6.1.1  Radiation  damage 

All  cells  in  the  lymphopoiesis  model  are  considered  radiosensitive,  and  once  damaged,  decay 
at  rate  fi. 

6.1.2  Feedback  loops 

The  only  feedback  loop  included  in  the  lymphopoiesis  model  is  the  generic  mediator  that 
stimulates  the  repopulation  rate  of  X\.  The  mediator’s  concentration  decreases  as  the  pop- 
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ulation  of  lymphopoietic  cells  increases.  Because  the  action  of  the  mediator  is  fast,  its 
concentration  can  be  written  as  a  function  of  cell  concentration,  and  thus  the  rate  of  re¬ 
population  (B)  can  also  be  written  as  a  function  of  cell  concentration  (see  Section  3.2  for 
derivation): 


1  +  f3(x2d  +  xf  +  62(x2d  +  xd)  +  03(x3d  +  £3)) 

The  sources  for  rates  between  compartments  (7,  5),  the  rates  of  cell  decay  (0,  /i),  ra¬ 
diosensitivity  parameters  (D®,  ni,  D °,  n2,  D3,  n3),  and  other  parameters  (a,  02,  03)  are  given 
in  Section  6.3. 

6.2  Mathematical  Model 

The  mathematical  model  of  lymphopoiesis  is  described  by  a  set  of  initial  conditions  and  dif¬ 
ferential  equations  that  describe  the  change  in  cell  concentrations  with  time.  The  description 
of  each  cell  type  is  given  in  Table  6.1. 

Table  6.1:  Cell  concentrations  in  each  compartment  of  the  lymphopoiesis  model 


Variable  name  Description 

x\d  Undamaged  cells  in  the  mitotic  compartment  X\ 

x2d  Undamaged  cells  in  the  post  mitotic  compartment  X2 

x3d  Undamaged  lymphocytes  in  the  blood  X3 

xd  Damaged  cells  in  the  mitotic  compartment  X] 

xd  Damaged  cells  in  the  post  mitotic  compartment  A2 

xd  Damaged  lymphocytes  in  the  blood  X3 

x™d  Concentration  of  cells  in  the  Ah  compartment  at  equilibrium  (i=l,2,3) 


At  equilibrium  cell  concentrations  in  the  Ad,  X2l  and  Ad  compartments  are  equal  to  x2d, 
x2d,  and  x3d,  respectively.  Following  acute  exposure,  cells  immediately  transition  into  the 
damaged  compartment  depending  on  their  radiosensitivity.  The  following  equations  describe 
how  cell  concentrations  change  following  an  acute  exposure  of  dose  D  and  determine  initial 
conditions  for  the  temporal  model. 


xf  =  xf  (l  -  (1  -  e~D/D°i)ni )  (i  =  1,  2,  3)  (72) 

xd  =  xf  (l-e~D/D^yi  (i  =  1,2,3)  (73) 

Once  radiation  exposure  occurs,  a  series  of  differential  equations  are  used  to  describe 
the  compartmental  transitions,  repopulation  and  decay  rates,  and  the  feedback  mechanisms 
shown  in  Figure  6.1.  The  following  equations  describe  the  model  where  Eq.  72-73  are  used 
to  calculate  initial  conditions.  Xi  represents  dxi/dt. 
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where 


xf  =  Bxf  -  "fxf 

(74) 

xf  =  7x-f  -  Sxf 

(75) 

xf  =  Sxf  -  *l)xf 

(76) 

xdi  =  (i  =  1,  2,  3) 

(77) 

1  +  /3(xf  +  xf  +  02(x%d  +  xf)  +  03(xf  +  xf)) 

6.3  Radiation  Model  Parameter  Sources 

Lymphopoiesis  parameters  were  determined  either  using  experimental  data  or  through  op¬ 
timization.  The  optimized  parameters  include  7,  #2,  and  D The  following  gives  a  detailed 
explanation  of  the  source  of  all  other  parameters  in  the  human  lymphopoiesis  model. 

6.3.1  Maximum  mitotic  repopulation  rate  (a) 

The  maximum  mitotic  repopulation  rate  was  set  to  2.0  d^1  which  results  in  a  minimum 
division  time  of 


_  ln( 2)  24 h 

~  2d-1  '  a 


8.3  h 


(79) 


In  adult  humans,  cells  can  divide  up  to  about  twice  a  day,  and  the  DNA  synthesis  phase 
only  takes  approximately  8  h  (Andreeff  et  al.  2000).  Thus,  a  max  repopulation  rate  of  2.0 
d_1  is  biologically  reasonable. 


6.3.2  Rate  of  post-mitotic  progenitor  maturation  (6) 

The  rate  of  post-mitotic  progenitor  maturation  was  determined  by  estimating  the  number 
of  post-mitotic  progenitors  in  the  bone  marrow  and  thymus  and  deriving  the  lymphocyte 
turnover  rate  by  this  estimated  post-mitotic  cell  concentration.  In  humans,  the  estimated 
percentage  of  lymphocytes  in  the  bone  marrow,  thymus,  and  blood  is  10.9%,  10.9%,  and 
2.2%,  respectively  (Trepel  1974).  Assuming  lymphocytes  in  the  bone  marrow  and  thymus 
comprise  X2,  this  leads  to  an  X3/X2  ratio  of  0.10.  At  equilibrium: 

S  =  0  —  (80) 

X2 

Thus,  5  is  equal  to  0.100.  In  the  next  section,  the  value  of  0  in  humans  is  determined  to  be 
0.01  d-1  meaning  5  equals  approximately  0.001  d_1  in  humans. 
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6.3.3  Rate  of  lymphocyte  clearance  from  blood  (ip) 

In  humans,  the  fractional  replacement  rates  of  CD4+  and  CD8+  T  cells  are  0.008  ±  0.005  d^1 
and  0.009  ±  0.013  d_1,  respectively  (Hcllerstein  et  al.  1999).  At  equilibrium  in  our  model,  the 
rate  of  replacement  is  equivalent  to  the  rate  of  decay.  B  cells  in  young  individuals  (age=19- 
34)  were  found  to  have  a  mean  proliferation  rate  of  0.0193  ±  0.0103  d-1  and  a  mean  clearance 
rate  of  0.0394  ±  0.0339  d-1  (Macallan  et  al.  2005).  Because  our  model  does  not  take  into 
account  lymphocyte  proliferation,  the  value  of  ip  for  B  cells  can  be  obtained  by  subtracting 
the  clearance  rate  from  the  proliferation  rate.  This  gives  a  decay  rate  applicable  to  our 
model  of  0.0146  d-1.  This  analysis  leads  to  an  approximate  decay  rate  for  all  lymphocytes 
of  0.01  d"1. 

6.3.4  Effect  of  damaged  states  on  mediator  levels  (/r) 

The  death  rate  of  damaged  cells  was  set  to  1.0  d-1.  In  Smirnova’s  original  model,  the 
moderately  and  heavily  damaged  compartments  had  death  rates  of  0.5  d-1  and  6.0  d^1  to 
simulate  mitotic  and  interphase  death,  respectively  (Smirnova  2011).  Here,  the  damaged 
compartment  contains  both  cells  that  die  during  mitosis  and  interphase.  Thus,  the  value  of 
1.0  d-1  provides  an  approximate  average  lifespan  of  1.0  d. 

6.3.5  Radiosensitivity  of  post-mitotic  precursors  ( D 2) 

The  radiosensitivity  of  post- mitotic  precursors  was  assumed  to  be  equal  to  that  of  mitotic  pre¬ 
cursors:  D®  =  D\.  In  Smirnova’s  murine  model  of  lymphopoiesis  the  mitotic  and  post- mitotic 
precursors  have  the  same  radiosensitivity  (Smirnova  2011),  and  we  use  that  assumption  here 
in  the  human  model. 

6.3.6  Radiosensitivity  of  lymphocytes  (D®,  n3) 

The  radiosensitivity  of  lymphocytes  in  humans  was  determined  by  fitting  the  survival  curve 
to  in  vitro  normalized  lymphocyte  survival  data  4  days  post  irradiation  (Kwan  and  Norman 
1977).  The  optimized  value  of  D®  was  6.22  Gy,  and  the  optimized  value  of  n3  was  0.40. 

6.3.7  Targets  per  progenitor  (n\,  n2) 

The  average  number  of  targets  per  cell  in  the  lymphopoiesis  model  was  assumed  to  be  1  for 
the  mitotic  (rq  =  1)  and  post-mitotic  progenitor  cells  (n2  =  1).  Using  this  assumption,  we 
are  simplifying  our  radiation  model  to  the  single-target,  single-hit  model  used  by  Smirnova 
(Smirnova  2011).  Using  this  model  of  radiation  damage,  we  are  able  to  accurately  match 
cell  survivability  following  radiation  exposure. 

6.3.8  Parameter  dependencies  (d3) 

The  amount  of  mediator  decay  caused  by  x3  and  x2  cells  is  assumed  to  be  identical  [02  =  #3). 
To  our  knowledge  there  is  no  experimental  data  available  to  determine  the  relationship 
between  02  and  #3.  Perturbations  to  these  parameters  affect  the  model  output  in  similar 
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Table  6.2:  Optimization  data  summary  for  lymphopoiesis  model 


Incident /Subject 

Reference 

Dose  (Gy) 

Baseline  (lO^L  x) 

Los  Alamos  1945-2 

Hempclmann  et  al.  1952 

0.12 

2.45“ 

UT  CARL  1971 

Andrews  et  al.  1961* 

2.6 

2.45“ 

Argonne  1952-1 

Bond  et  al.  1965 

1.5 

2.45“ 

Y-12  1958-A 

Bond  et  al.  1965 

4.02 

2.45“ 

Y-12  1958-B 

Andrews  et  al.  1961* 

2.97 

2.45“ 

Y-12  1958-C 

Andrews  et  al.  1961* 

3.73 

2.45“ 

Y-12  1958-E 

Gusev  et  al.  2010* 

2.59 

2.45“ 

Vinca  1958-B 

Jammet  et  al.  1959 

2.97 

2.45“ 

Vinca  1958-H 

Jammet  et  al.  1959 

4.37 

2.45“ 

Vinca  1958-M 

Jammet  et  al.  1959 

5.74 

2.45“ 

Vinca  1958-G 

Jammet  et  al.  1959 

5.35 

2.45“ 

Vinca  1958-D 

Jammet  et  al.  1959 

5.40 

2.45“ 

Vinca  1958-V 

Jammet  et  al.  1959 

6.11 

2.45“ 

Cancer  patients  1958 
-mean  of  30 

Miller  et  al.  1958 

1.0 

1.395b 

Cancer  patients  1958 
-mean  of  12 

Miller  et  al.  1958 

1.5 

1.3856 

Cancer  patients  1958 
-mean  of  18 

Miller  et  al.  1958 

2.0 

1.35b 

Yarmonenko  1988-Z 

Yarmonenko  1988 

9.8 

2.45“ 

Lockport  1960-A 

Bond  et  al.  1965 

3.0 

2.45“ 

Lockport  1960-B 

Bond  et  al.  1965 

3.0 

2.45“ 

“Mean  platelet  concentration  observed  in  humans  (Valentin  2002) 
bMean  of  pre-radiation  data  points 
*Data  provided  by  Ron  Goans 

Additional  details  on  the  human  case  studies  provided  in  Appendix  A. 

ways.  Thus,  when  both  are  included  in  the  optimization,  the  model  is  not  identifiable.  Here, 
the  model  was  simplified  by  setting  these  parameters  equal  to  each  other. 

6.4  Optimization  Data 

For  the  radiation  model  optimization,  data  was  collected  from  multiple  sources.  Table  6.2 
provides  a  summary  of  the  radiation  data  used  in  optimizing  the  human  lymphopoiesis 
model.  Information  on  the  data  reference,  dose,  and  baseline  values  is  given  for  all  case 
studies.  Data  from  19  case  studies  were  used  in  the  optimization,  12  of  which  contained 
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data  points  up  to  and  beyond  30  days  post-radiation  exposure.  The  subjects  received  doses 
ranging  from  0.12  to  9.8  Gy. 

6.5  Results:  Optimization  and  MCMC  Analysis 

Table  6.3  gives  a  complete  list  of  parameter  values,  descriptions,  and  bases.  The  optimized 
parameter  values  and  corresponding  95%  confidence  intervals  are  given  in  Table  6.4.  Also 
provided  are  the  initialization  value,  the  lower  bound,  and  the  upper  bound  placed  on  each 
parameter  prior  to  running  the  optimization.  Figure  6.2  shows  the  optimized  model  com¬ 
pared  to  select  data  used  for  the  optimization.  Each  dataset  comes  from  the  specified  case 
study/patient.  The  remainder  of  the  optimization  data  is  compared  to  the  model  output 
either  in  the  Discussion  (Section  7;  Figure  7.5b)  or  in  Appendix  B.l.  For  details  on  the 
MCMC  results  and  identihability  analysis  see  Appendix  C. 


Table  6.3:  Biological  parameter  values  for  lymphopoiesis 


Parameter 

Description 

Value 

Section  (References) 

a 

Maximum  rate  of  mitotic  pro¬ 
genitor  cell  repopulation 

2  d-1 

6.3.1  (Andreeff  et  al.  2000) 

7 

Rate  of  mitotic  progenitor  cell 
maturation 

0.06  d”1 

Optimized 

5 

Rate  of  post-mitotic  progenitor 
cell  maturation 

0.001  d”1 

6.3.2  (Trepel  1974) 

Rate  of  lymphocyte  clearance 

0.01  d"1 

6.3.3  (Hellcrstein  et  al.  1999; 

from  the  blood 

Macallan  et  al.  2005) 

h 

Rate  of  damaged  cell  death 

1  d”1 

6.3.4 

D°i 

Determines  fraction  of  damaged 
X\  cells 

3.98  Gy 

Optimized 

D° 

Determines  fraction  of  damaged 
x2  cells 

3.98  Gy 

6.3.5  (Smirnova  2011) 

n° 

u?> 

Determines  fraction  of  damaged 
x3  cells 

6.22  Gy 

6.3.6  (Kwan  and  Norman  1977) 

nt 

Number  of  targets  per  cell  X; 
(*  =  1,2) 

1 

6.3.7 

n3 

Number  of  targets  per  x3  cell 

0.40 

6.3.6  (Kwan  and  Norman  1977) 

e2 

Decay  rate  of  mediator  due  to 
x2  cells  relative  to  decay  rate 
due  to  X\  cells 

2.67 

Optimized 

o3 

Decay  rate  of  mediator  due  to 
x3  cells  relative  to  decay  rate 
due  to  x\  cells 

2.67 

6.3.8 
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Table  6.4:  Lymphopoiesis  optimized  parameters  and  confidence  intervals 


Parameter 

Optimized  Value  2.5%“ 

97. 5%6 

Initial  Value  (Lower, Upper)c 

7  (d-1) 

0.0605  0.0534 

0.1243 

0.0447  (0.001,  1.99) 

e2 

2.6705  0.5312 

3.167 

19.9875  (0.001,  20) 

D°i  (Gy) 

3.9769  2.8992 

4.4005 

4.1128  (0.1,  20) 

“Lower  end  of  parameter  95%  confidence  interval 
&Upper  end  of  parameter  95%  confidence  interval 

initialization  value,  lower  (min),  and  upper  (max)  bound  given  to  parameter  during  optimization 


6.6  Validation 

To  validate  the  human  lymphopoiesis  model,  the  model  outputs  were  compared  to  data  on 
lymphocyte  count  following  radiation  in  humans  that  were  not  included  in  the  optimization 
(Hempelmann  et  ah  1952;  Bond  et  ah  1965;  Hirama  et  ah  2003;  Stavern  et  ah  1985;  Mettler 
et  ah  2001;  Mettler  2001).  Select  validation  results  are  shown  in  Figure  6.3.  An  exception 
to  this  is  the  Mol  1965  case  study  which  was  a  highly  non-uniform  exposure.  The  Mol 
exposure  is  included  in  the  validation  figure  to  provide  a  comparison  with  the  Sarov  1963 
case  study,  which  is  intended  for  validation,  who  received  the  same  estimated  dose  of  5.5  Gy. 
The  remainder  of  the  validation  results  are  given  in  Appendix  B.2. 
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CO 


(a)  Dose:  1.0  Gy 


(b)  Dose:  2.97  Gy  ±  25% 


CO 


CO 


(c)  Dose:  4.37  Gy  ±  15%  (d)  Dose:  9.8  Gy 

Figure  6.2:  Lymphopoiesis  model  compared  to  select  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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CO 


(a)  Dose:  0.65  Gy 


CO 


(c)  Dose:  5.5  Gy 


CO 


(b)  Dose:  3.6  Gy 


CO 


Time  following  radiation  (days) 

(d)  Dose:  22.5  Gy 


Figure  6.3:  Lymphopoiesis  model  compared  to  select  validation  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line.  Output  initialized  at 
the  upper  and  lower  end  of  the  lymphocyte  range  is  shown  by  dotted  lines. 
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7  Discussion 


Three  hematopoietic  mechanistic  models  which  simulate  blood  cell  kinetics  following  acute 
radiation  exposure  in  humans  are  presented.  These  models  have  been  optimized  with  and 
validated  against  available  case  study  data.  The  aim  of  our  modeling  work  was  to  develop 
models  that  represent  the  general  hematopoietic  response  in  humans  after  radiation  exposure 
without  treatment.  Our  models  predict  the  average  response  in  humans,  provide  realistic 
estimates  of  blood  cell  counts,  and  capture  the  important  trends  observed  in  blood  cell 
dynamics  following  radiation  exposure  at  a  range  of  dose  levels. 

Several  challenges  arise  when  optimizing  and  validating  the  models.  For  both  these 
procedures,  we  compare  human  data  to  model  simulations.  However,  a  number  of  factors 
impact  the  relevancy  and  reliability  of  these  comparisons: 

1.  Baseline  blood  cell  counts  are  not  available  in  the  majority  of  case  studies.  For  each 
blood  cell  type,  there  is  a  large  range  for  what  this  baseline  value  could  be. 

2.  Uncertainty  in  acute  radiation  doses  due  to: 

(a)  Variable  reporting  of  free-in-air  versus  midline  tissue  doses 

(b)  Uncertainty  in  physical  dose  reconstruction  resulting  in  a  large  dose  range 

(c)  Uncertainty  in  the  dose  extrapolation  for  different  radiation  types  (i.e.  neutrons 
and  x-rays) 

(d)  Non-uniformity  of  dose  delivered 

(e)  Dose  rate  and  duration;  our  models  assume  that  the  dose  is  received  instanta¬ 
neously 

3.  Treatments,  such  as  transfusion  and  cytokine  therapy,  may  impact  the  observed  blood 
cell  kinetics. 

4.  Inter- individual  variabilities  are  observed  in  blood  cell  kinetics  and  cell  radiosensitivi¬ 
ties. 

To  perform  comparisons  to  the  dimensionless  model  output,  data  is  divided  by  a  base¬ 
line  value  specific  to  each  individual.  Healthy  humans  have  a  large  range  of  possible  blood 
cell  concentrations;  thus,  an  accurate  baseline  value  is  essential  for  comparisons.  For  opti¬ 
mization  if  the  baseline  value  is  incorrect,  the  calculated  residual  values  may  be  artificially 
large  or  small.  For  granulocytes  and  lymphocytes  when  baseline  cell  counts  specific  to  the 
individual  were  not  available,  the  average  concentration  observed  in  humans  was  used  as  the 
baseline  (Valentin  2002).  For  platelets,  because  there  is  a  delayed  response  following  radi¬ 
ation  exposure,  platelet  counts  obtained  within  three  days  of  radiation  exposure  were  used 
for  determining  baseline  values.  If  values  before  three  days  were  not  available,  the  average 
concentration  observed  in  humans  was  used  (Valentin  2002). 

To  address  concerns  about  the  uncertainty  in  radiation  dose  estimates,  case  studies  were 
reviewed  in  detail.  An  effort  was  made  to  identify  studies  that  involved  acute  exposures  that 
were  primarily  uniform  with  reliable  dose  estimates.  Nevertheless,  all  exposures  result  in 
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some  non- uniformity  and  uncertainty  in  the  reported  doses.  Furthermore,  due  to  limitations 
in  available  data,  some  comparisons  are  made  in  spite  of  known  issues  with  the  data.  In 
these  cases,  the  conditions  must  be  taken  into  consideration  when  comparing  with  model 
outputs. 

An  effort  was  also  made  to  identify  treatment  protocols  for  each  case  study.  While  most 
patients  did  receive  some  treatment,  our  main  concern  was  to  identify  blood  transfusions, 
stem  cell  transplants,  and  cytokine  therapy  since  these  procedures  can  have  a  significant 
impact  on  the  observed  blood  cell  kinetics.  Again,  due  to  limitations  in  available  data,  some 
comparisons  are  performed  with  data  from  treated  subjects  or  from  subjects  where  treatment 
information  is  unknown.  The  use  of  antibiotics  for  infections  was  not  controlled  for  in  the 
data.  These  factors  may  impact  blood  cell  kinetics,  and  we  plan  to  address  this  question 
in  future  work  by  analyzing  experimental  animal  data  in  which  other  conditions  are  more 
tightly  controlled. 

Individual  variability  in  cell  kinetics  and  radiosensitivity  also  impacts  comparisons  of  the 
model  to  individual  data.  Parameter  values  obtained  from  the  literature  are  often  presented 
as  a  mean  value  with  an  error  measurement  that  provides  information  on  the  parameter 
variance.  For  example  the  granulocyte  transit  time  through  the  mitotic  compartment  is 
6.4 ±0.3  d  (Price  et  al.  1996).  Also,  when  estimating  parameter  values  through  optimization, 
variance  estimates  on  the  optimized  parameter  values  were  obtained.  However,  in  the  current 
models,  only  the  mean  parameter  value  is  used  and  information  on  the  error  associated  with 
the  parameter  value  is  lost.  Using  a  population  based  modeling  approach  which  performs 
simulations  accounting  for  these  parameter  variations  may  help  understand  and  quantify  the 
variance  in  output.  In  the  future,  we  plan  to  conduct  population  analyses  and  Monte  Carlo 
simulations  to  predict  the  range  and  likelihood  of  model  outputs. 

7.1  Thrombopoiesis 

The  thrombopoiesis  model  accurately  captures  important  dose-dependent  trends  observed 
following  acute  radiation  exposure.  Figure  7.1a  shows  the  thrombopoiesis  model  simulations 
at  increasing  radiation  dose  levels,  and  Figure  7.1b  compares  simulations  to  observational 
data.  The  purpose  of  these  comparisons  is  to  verify  that  trends  observed  in  the  data  are 
reflected  by  the  model. 

At  low  acute  radiation  doses,  there  is  a  plateau  phase  prior  to  a  decrease  in  platelet  count. 
The  duration  of  this  phase  represents  the  time  needed  for  megakaryocyte  maturation.  As  the 
radiation  dose  increases,  megakaryocytes  start  becoming  damaged  and  the  plateau  phase  is 
eliminated.  The  platelet  counts  reach  a  nadir  between  20  and  30  days  following  exposure  and 
subsequently  begin  to  recover.  This  recovery  results  in  an  overshoot  where  platelet  levels 
temporarily  rise  above  normal  values.  The  significant  trends  observed  in  the  data  as  the 
radiation  dose  increases  are: 

1.  Platelet  counts  show  a  more  drastic  initial  decline,  and  the  plateau  phase  is  eliminated. 

2.  The  nadir  is  lower  but  occurs  roughly  around  the  same  time. 

The  model  of  thrombopoiesis  is  able  to  capture  these  trends. 
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(a) 


(b) 


Figure  7.1:  Thrombopoiesis  simulations 


One  challenge  encountered  in  optimization  of  the  thrombopoiesis  model  was  reliably 
capturing  the  depth  of  the  nadir.  However,  although  there  is  a  clear  correlation  between  the 
nadir  depth  and  the  radiation  dose;  in  some  case  study  pairs,  the  depth  of  the  nadir  is  not 
lower  in  the  case  study  with  the  larger  radiation  dose.  For  example,  as  shown  in  Figure  7.2a, 
the  minimum  normalized  platelet  count  for  the  2.6  Gy  case  study  from  UT  CARL  (Goans 
2012)  was  lower  than  that  observed  in  the  3.51  Sv  case  study  B  from  the  Vinca  incident 
(Jammet  et  al.  1959;  Pesic  2012).  This  inconsistency  could  arise  from  individual  variability 
in  the  response  to  radiation  or  result  from  uncertainty  in  the  reported  dose.  Similarly,  in 
the  5.35  Gy  case  study  G  from  Vinca  (Jammet  et  al.  1959;  Pesic  2012),  platelet  counts 
decline  less  than  in  the  Y-12  case  A  subject  (Bond  et  al.  1965;  Mettler  et  al.  2001;  Kerr 
and  Tankersley  2006)  who  received  an  estimated  4.02  Gy.  A  third  comparison,  shown  in 
Figure  7.2b,  shows  the  Vinca  V  subject  who  received  an  estimated  6.11  Gy  dose  and  Subject 
Z  from  Yarmonenko  1988  who  received  an  estimated  9.8  Gy  dose.  The  Vinca  V  subject 
has  a  more  drastic  decline  in  platelet  counts  immediately  following  radiation  exposure  when 
compared  with  Subject  Z;  however,  we  would  expect  to  see  a  faster  decline  in  subject  Z  since 
that  subject  received  a  higher  dose.  Although  all  these  comparisons  are  performed  using  the 
normalized  data,  the  same  issues  arise  when  performing  the  comparisons  using  the  absolute 
platelet  counts.  Such  discrepancies  in  the  case  study  data  help  explain  why  some  of  the 
simulation  comparisons  vary  significantly. 

Validation  blood  cell  kinetic  data  following  radiation  exposure,  presented  in  Figure  4.4 
and  Figure  B.4,  compare  favorably  with  the  simulation  results.  Considerable  validity  of  the 
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Figure  7.2:  Thrombopoiesis  case  study  comparisons 


thrombopoiesis  model  is  also  gained  through  the  thrombopheresis  data  comparisons  (Figure 
4.5).  In  four  separate  studies,  two  of  which  had  mean  data  from  multiple  subjects,  the 
model  was  able  to  accurately  predict  the  platelet  count  recovery  following  platelet  removal. 
These  comparisons  reassure  us  that  the  kinetic  parameters  obtained  from  the  literature  and 
through  optimization  are  accurate.  Additionally,  using  model  parameters  to  calculate  the 
relative  size  of  the  MK  compartment  compared  to  the  platelet  compartment,  as  described 
in  Section  4.6.1,  led  to  relative  compartment  sizes  that  are  in  agreement  with  experimental 
data  (Levine  1980;  Valentin  2002;  Starr  and  Taggart  1989). 

While  large  variances  are  observed  in  some  of  the  case  study  data,  a  general  dose  re¬ 
sponse  has  been  established  that  reliably  predicts  the  trends  observed  in  platelet  decline 
after  radiation  exposure.  Further  validation  for  the  structure  and  parameter  values  of  the 
model  is  added  through  both  the  thrombopheresis  data  comparisons  and  the  compartment 
ratio  analysis. 


7.2  Granulopoiesis 

The  granulopoiesis  model  accurately  predicts  observed  trends  in  circulating  granulocyte 
count  following  acute  radiation  exposure.  Figure  7.3a  shows  the  granulopoiesis  model  simu¬ 
lations  at  increasing  radiation  dose  levels,  and  Figure  7.3b  compares  similar  simulations  to 
observational  data.  The  purpose  of  these  comparisons  is  to  verify  that  trends  observed  in 
the  data  are  reflected  by  the  model. 
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Figure  7.3:  Granulopoiesis  simulations 


Following  radiation  exposure  there  is  an  initial  decrease  in  granulocyte  count  until  the 
first  nadir  is  reached.  The  system  then  appears  to  begin  recovery;  this  phenomena  is  known 
to  be  an  abortive  rise  in  granulocyte  count.  At  between  12  and  20  days,  the  granulocyte 
count  decreases  again  until  a  second  nadir  is  reached.  The  significant  trends  observed  as  the 
radiation  dose  increases  are: 

1.  The  time  of  abortive  rise  occurs  earlier. 

2.  The  nadirs  before  and  after  the  time  of  abortive  rise  are  lower. 

3.  The  recovery  of  the  system  occurs  earlier. 

The  model  of  granulopoiesis  is  able  to  predict  these  dose-dependent  trends. 

As  seen  in  the  thrombopoiesis  work,  in  some  case  studies  individual  variability  or  envi¬ 
ronmental  factors  clearly  impact  general  trends,  making  comparisons  difficult.  For  example, 
the  subject  in  the  UT  CARL  case  study  (Goans  2012)  experienced  a  lower  radiation  dose 
of  2.6  Gy  compared  with  the  3.51  Gy  dose  received  by  subject  B  from  the  Vinca  incident 
(Jarnmet  et  al.  1959;  Pesic  2012).  Figure  7.4a  shows  a  comparison  of  the  two  datasets.  As 
expected  the  first  granulocyte  nadir  around  day  10  is  higher  in  LIT  CARL  subject  than  in  the 
Vinca  B  subject.  However,  for  the  second  nadir,  the  UT  CARL  subject  has  a  lower  minimum 
granulocyte  concentration  compared  to  the  Vinca  B  subject.  If  this  were  a  trend  observed 
between  other  datasets,  we  would  consider  revising  the  model  structure.  However,  because 
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Figure  7.4:  Granulopoiesis  case  study  comparisons 


this  switch  in  minimum  nadir  occurs  rarely,  we  attribute  it  to  the  individual  variability  and 
treatment  environment  of  the  subjects  rather  than  to  an  error  in  our  model  structure. 

Figure  7.4b  shows  a  second  comparison  between  subjects  G  and  D  from  the  Vinca.  Both 
subjects  received  a  radiation  dose  of  approximately  5.4  Gy.  The  granulocyte  kinetics  of  sub¬ 
ject  G  would  suggest  a  lower  dose  compared  with  subject  D.  Since  both  these  subjects  have 
a  radiation  dose  uncertainty  of  ±15%,  we  attribute  the  differences  in  granulocyte  kinetics  to 
uncertainty  in  the  dose.  However,  we  can  not  say  for  sure  that  subject  D  received  a  higher 
dose  because  treatment  and  individual  variations  may  still  play  a  role. 

In  many  subjects,  an  early  abortive  rise  in  granulocyte  counts  immediately  following  ra¬ 
diation  is  observed.  The  current  model  does  not  predict  this  early  rise  which  is  predicted 
to  result  from  an  immediate  release  of  granulocytes  from  the  marginated  pool  (granulocytes 
adhering  to  the  walls  of  the  endothelium)  into  circulation.  The  X:i  compartment  of  the  gran¬ 
ulopoiesis  model  includes  both  circulating  and  marginated  granulocytes,  and  therefore  the 
model  does  not  depict  the  early  increase  in  granulocyte  count.  Because  the  feedback  mecha¬ 
nisms  in  the  model  depend  on  the  total  number  of  granulocytes,  the  release  of  granulocytes 
from  the  marginated  pool  does  not  affect  these  mechanisms.  When  making  comparisons  be¬ 
tween  the  model  and  circulating  granulocyte  data,  we  assume  that  the  ratio  of  marginated 
to  circulating  granulocytes  remains  constant.  Although  this  assumption  is  not  valid  immedi¬ 
ately  following  radiation  exposure,  once  the  ratio  of  marginated  to  circulating  granulocytes 
returns  to  normal,  the  current  model  is  still  valid.  Model  adjustments  to  account  for  this 
early  change  in  the  marginated  to  circulating  granulocyte  ratio  can  be  added  if  this  infor- 
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mation  is  needed  in  the  future. 

Factors  such  as  treatment,  non-uniform  exposures,  and  internal  contamination  limit  the 
robustness  of  our  validation  comparisons  for  the  granulocyte  model.  Since  the  granulocyte 
cell  line  responds  quickly  to  radiation,  assigning  accurate  baseline  values  for  the  case  studies 
is  difficult.  Therefore,  comparisons  are  presented  with  simulation  outputs  initialized  at  the 
population  mean  granulocyte  concentration  together  with  the  upper  and  lower  bounds  of 
the  granulocyte  concentration  observed  in  healthy  humans. 

Validation  data  presented  in  Figure  5.3  and  Figure  B.5  for  the  most  part  compare  fa¬ 
vorably  with  the  simulation  results,  with  a  few  exceptions.  Los  Alamos  subject  6  appears 
to  have  very  high  granulocyte  blood  cell  counts  during  the  entire  duration  of  the  study, 
suggesting  the  baseline  value  used  was  low  for  that  individual. 

7.3  Lymphopoiesis 

The  lymphopoiesis  model  accurately  captures  observed  trends  in  lymphocyte  count  follow¬ 
ing  acute  radiation  exposure.  Figure  7.5a  shows  the  lymphopoiesis  model  simulations  at 
increasing  radiation  dose  levels,  and  Figure  7.5b  compares  four  simulations  to  observational 
data.  The  purpose  of  these  comparisons  is  to  verify  that  trends  observed  in  the  data  are 
reflected  by  the  model. 
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Radiation  leads  to  a  rapid  decline  in  lymphocyte  counts  followed  by  a  slow  recovery.  The 
main  dose-dependent  trend  in  the  lymphopoiesis  model  is  the  rate  of  decline  increases  as 
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the  radiation  dose  becomes  greater.  Our  model  of  lymphopoiesis  is  able  to  accurately  depict 
how  the  lymphocyte  nadir  changes  as  the  dose  increases  and  the  subsequent  recovery. 

Once  again,  discrepancies  between  case  studies  and  general  dose-dependent  trends  exist; 
however,  these  issues  are  not  as  prominent  as  in  the  platelet  and  granulocyte  cases.  An 
example  discrepancy  is  shown  in  Figure  7.6.  Here,  the  Argonne  I  subject  who  received 
a  1.5  Gy  radiation  dose  is  compared  to  the  Y-12  1958  E  subject  who  received  a  2.59  Gy 
radiation  dose.  In  the  2.59  Gy  case  study,  although  the  normalized  values  for  the  lymphocyte 
concentration  reach  a  lower  value  sooner  than  in  the  Argonne  1.5  Gy  case,  at  later  time  points 
the  normalized  lymphocyte  concentration  from  the  2.29  Gy  case  study  is  either  equal  to  or 
greater  than  the  lymphocyte  concentration  from  the  Argonne  1.5  Gy  case.  The  lymphopoiesis 
model  output  matches  the  Argonne  1.5  Gy  study  better  than  the  Y-12  E  study,  suggesting 
the  issue  is  due  to  dose  uncertainty  or  individual  variability  in  the  Y-12  E  study. 
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Figure  7.6:  Lymphopoiesis  case  study  comparisons 


Validation  data  presented  in  Figure  6.3  and  Figure  B.6  compare  favorably  with  the  sim¬ 
ulation  results,  with  a  few  exceptions.  At  the  higher  doses,  the  lymphocyte  count  appears 
to  drop  faster  than  the  model  predicts  (Case  Study:  Tokai-mura  1999-B,  Shanghai-S,  and 
Tokai-mura  1999-A).  If  this  early  response  needs  to  be  captured,  future  work  would  involve 
including  the  rate  of  damaged  cell  death  in  the  optimization  procedure. 


7.4  Comparison  with  Expert  Predictions 

METREPOL  (Medical  Treatment  Protocols  for  Radiation  Accident  Victims  as  a  Basis  for  a 
Computerised  Guidance  System)  provides  an  assessment  of  radiation  injury  severity  based 
upon  the  hematopoietic  response  (Fliedner,  Friesecke,  et  al.  2001).  This  severity  ranking 


system  was  developed  by  an  international  group  of  selected  experts  and  is  based  on  both 
experimental  data  and  the  analysis  of  data  from  radiation  victims.  There  are  four  possible 
severity  grades  (HI:  Mild  damage;  H2:  Moderate  damage;  H3:  Severe  damage;  H4:  Fatal 
damage)  which  are  not  based  on  radiation  dose  but  on  the  ensuing  symptoms. 

To  verify  that  our  hematopoietic  models  are  in  agreement  with  the  predicted  hematopoi¬ 
etic  responses  from  METREPOL  and  to  correlate  these  responses  with  a  dose,  we  developed 
overlays  of  our  model  on  the  METREPOL  predictions  for  the  H2  and  H3  severity  levels  (Fig¬ 
ure  7.7).  For  an  accurate  comparison,  our  hematopoietic  model  outputs  are  shown  for  the 
normal  range  of  blood  cell  counts  as  given  in  the  METREPOL  report  (Fliedner,  Friesecke, 
et  al.  2001).  Interestingly,  for  the  same  severity  ranking  the  hematopoietic  models  were 
initialized  at  different  acute  radiation  doses  to  compare  to  H2  and  H3  METREPOL  pre¬ 
dictions.  Since  hematopoietic  cell  lines  respond  differently  to  radiation,  the  H1-H4  severity 
levels  listed  in  METREPOL  may  correspond  to  different  doses  for  the  same  severity  level. 
For  example,  to  match  the  predicted  H2  response,  the  lymphopoiesis  model  was  run  with  a 
radiation  dose  of  3  Gy  while  the  granulopoiesis  model  was  run  at  7  Gy. 

In  general,  the  model  predictions  match  the  trends  given  by  METREPOL.  The  platelet 
recovery  predicted  by  our  model  is  slightly  faster  and  leads  to  an  overshoot  not  predicted  in 
the  METREPOL  hematopoietic  response;  however,  this  overshoot  is  observed  in  several  of 
the  individual  case  studies  examined  in  this  report.  As  examples  see  case  studies  Argonne 
1952-1,  Y-12  1958-A,  and  Yarnronenko  1988-Z  shown  in  Figure  4.3.  The  granulopoiesis  model 
has  an  abortive  rise  at  a  similar  time  as  the  METREPOL  predictions,  however  at  higher 
severity  our  granulopoiesis  model  predicts  an  earlier  recovery  due  to  the  shortened  time  of 
abortive  rise.  METREPOL  predicts  a  later  granulocyte  recovery  at  higher  doses.  Based  on 
the  dose-response  data  shown  in  Figure  7.3b,  an  early  recovery  of  granulocytes  at  higher 
doses  is  not  unfeasible. 

In  a  second  set  of  expert  predictions,  Annex  G  of  the  UNSCEAR  1988  report  presents 
a  schematic  picture  of  the  average  time  course  of  blood  cell  counts  deduced  from  accidental 
human  exposures  (UNSCEAR  1988).  Figure  7.8  shows  a  comparison  of  our  hematopoietic 
models’  outputs  with  this  expert  prediction.  The  UNSCEAR  report  presents  the  predicted 
blood  cell  responses  at  doses  ranging  from  0  to  20  Gy  divided  into  four  groups.  Therefore, 
we  plot  our  hematopoietic  model  outputs  as  a  shaded  region  based  on  each  dose  range  and 
overlay  the  predicted  UNSCEAR  response.  Although  the  models  are  not  in  exact  agreement 
with  the  UNSCEAR  predictions,  the  general  trends  are  supported  and  the  dose  response 
matches  fairly  well. 

7.5  Treatment  Effects 

The  purpose  of  our  models  is  to  show  hematopoietic  dynamics  following  untreated  radiation 
exposure;  however,  in  many  case  studies  treatments  were  administered  following  acute  radi¬ 
ation  exposure.  Comparing  treated  case  studies  to  the  model  prediction  provides  hypotheses 
for  how  treatments  affect  the  hematopoietic  response.  Here,  we  discuss  one  case  study  from 
Soreq  1990  (IAEA  1993;  Mettlcr  2001)  while  other  case  study  examples  involving  treatment 
are  presented  in  Appendix  B.2.  The  subject  from  Soreq  received  a  whole-body  exposure  of 
15  ±  5  Gy.  He  was  treated  with  G-CSF  one  day  following  exposure,  and  platelet  transfu¬ 
sions  were  given  when  the  platelet  count  was  less  than  20  •  10 3/iL”1.  Platelet  counts  first 
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Figure  7.7:  Comparison  of  hematopoietic  model  outputs  to  METREPOL 
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Figure  7.8:  Comparison  of  hematopoietic  model  outputs  (shaded  regions  for 
each  dose  range)  with  expert  opinions  on  blood  cell  dynamics  following  radiation 
exposure  (solid  and  dashed  lines;  UNSCEAR  1988) 
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Figure  7.9:  Hematopoietic  dynamics  following  15  Gy  (10-20  Gy)  and  treatment. 
Model  output  at  the  specified  dose  (range)  is  delineated  by  a  black  line  (shaded 
region).  If  a  generic  baseline  value  was  used,  output  initialized  at  the  upper  and 
lower  end  of  the  healthy  blood  cell  range  is  shown  by  dotted  lines. 
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decreased  below  this  point  on  the  9th  day  following  radiation  exposure.  Figure  7.9  shows 
the  kinetic  blood  cell  count  data  from  the  Soreq  subject.  The  granulocyte  data  suggests 
that  G-CSF  treatment  causes  a  rapid  recovery  of  the  granulocytes,  much  faster  than  what 
the  model  predicts.  The  peripheral  platelet  and  lymphocyte  concentrations  show  temporary 
rises,  which  might  be  due  to  treatment  as  well.  However,  it  would  take  more  data  collection 
and  analysis  to  verify  this. 

7.6  Partial  Body  Effects 

Several  case  studies  were  obtained  where  very  non-uniform  or  partial  body  exposures  oc¬ 
curred.  These  include  Los  Alamos  1945  subjects  1  and  3  (Hempelmann  et  al.  1952;  Hernpel- 
mann  1961;  Hoffman  and  Hempelmann  1957;  Mettler  et  al.  2001)  and  subject  K  from  Yarmo- 
nenko  1988.  Here  we  present  one  case  study,  Los  Alamos  subject  1,  who  received  an  acute 
localized  radiation  exposure.  Other  case  study  examples  involving  non-uniform  exposures 
are  presented  in  Appendix  B.2.  For  Los  Alamos  subject  1,  the  radiation  dose  was  estimated 
as  12.75  Sv  with  a  range  from  5.5-20  Sv.  Figure  7.10  shows  a  comparison  between  the  data 
and  model  output  for  all  three  cell  lines.  In  each  case,  the  blood  cell  kinetics  appear  to 
be  less  severe  than  what  the  models  predict.  We  would  expect  this  since  in  partial  body 
exposure  some  of  the  bone  marrow  is  spared.  HSCs  in  the  spared  bone  marrow  can  continue 
to  divide,  contribute  to  recovery,  and  help  alleviate  the  hematopoietic  radiation  effects. 

7.7  Current  and  Future  Effort 

While  the  hematopoietic  models  for  radiation  alone  provide  useful  insight  into  the  clinical 
course  of  blood  cell  kinetics,  mechanistic  models  also  provide  a  way  to  evaluate  the  impact  of 
additional  insults  on  hematopoiesis.  Therefore,  model  parameters  for  the  impact  of  burn  on 
hematopoiesis  are  also  under  development.  The  parameters  will  be  integrated  with  radiation 
models  to  simulate  the  effects  of  combined  radiation  and  burn  injury  in  humans.  These 
models  will  help  to  predict  the  synergistic  effects  of  combined  injury.  The  model  outputs 
will  also  be  used  to  facilitate  predictions  of  time  to  and  likelihood  of  mortality.  This  will  be 
done  by  developing  correlations  between  model  outputs  with  mortality.  Potential  outputs 
worth  examining  include  the  time  of  the  nadir  and  the  duration  in  which  blood  cell  counts 
are  below  specific  threshold  levels. 

The  hematopoietic  models  for  radiation  effects  could  also  be  used  in  the  future  to  explore 
the  overall  systemic  impact  of  partial  body  exposures  and  the  potential  decreased  mortality 
risk  afforded  by  partial  shielding  of  progenitor  (or  stem)  cells.  In  particular,  these  mechanistic 
models  can  help  sort  out  the  impact  of  global  mediator  response  on  blood  cell  kinetics  and 
the  associated  ability  of  the  hematopoietic  system  to  recover  after  inhomogeneous  radiation 
injury.  The  wealth  of  knowledge  gleaned  from  more  highly  controlled  animal  experiments 
will  be  useful  in  establishing  firm  parameters  for  a  partial  body  effects  model,  and  partial 
body  human  exposures  can  provide  data  for  optimization  and  validation  of  such  models. 

The  effect  of  treatment  on  hematopoietic  responses  after  radiation  exposure  is  also  an  area 
where  computational  approaches  may  be  useful.  In  separate  work,  we  have  integrated  the 
effect  of  the  cytokine  IL-12  on  hematopoietic  kinetics  on  these  models  based  on  key  changes 
observed  in  clinical  data  of  patients  that  received  IL-12  under  normal  circumstances  (i.e. 
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Figure  7.10:  Hematopoietic  dynamics  following  partial  body  exposure  of  12.75 
Sv  (5.5-20  Sv).  Model  output  at  the  specified  dose  (range)  is  delineated  by  a 
black  line  (shaded  region).  If  a  generic  baseline  value  was  used,  output  initialized 
at  the  upper  and  lower  end  of  the  healthy  blood  cell  range  is  shown  by  dotted 
lines. 
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without  radiation).  Since  G-CSF  is  currently  the  only  Food  and  Drug  Administration  (FDA) 
approved  therapeutic  for  the  hematopoietic  syndrome  associated  with  radiation  exposure, 
integration  of  its  effect  would  be  a  logical  next  step  in  our  modeling.  A  significant  amount  of 
data  on  G-CSF  is  available,  including  case  studies  of  its  use  with  highly  irradiated  persons, 
that  would  greatly  facilitate  parameterization  of  the  treatment  model.  Treatment  models 
can  aid  in  scenario  analyses  involving  large  numbers  of  casualties  and  limited  resources. 

Finally,  another  important  aspect  to  include  in  future  modeling  is  population  variability. 
In  our  current  work,  we  have  included  the  ability  to  adjust  baseline  blood  cell  counts  so  that 
the  effects  of  radiation  exposure  in  subpopulations  with  known  variances  in  baseline  cell 
counts  can  be  examined.  However,  a  more  thorough  population  modeling  approach  could  be 
taken  to  understand  the  population  variances  in  other  parameters  that  range  from  average 
blood  volumes  to  cell  turnover  and  regeneration  rates.  Integrating  known  parameters  for 
vulnerable  subpopulations  such  as  young  children  and  elderly  would  be  particularly  useful. 
This  type  of  approach  enables  model  users  to  understand  the  expected  variability  among  a 
large,  diverse  population  and  obtain  a  more  realistic  picture  of  the  health  effects  expected 
after  a  large  scale  nuclear  or  radiological  incident. 
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8  Conclusion 


Hematopoietic  models  describing  the  effects  of  radiation  on  thrombopoiesis,  granulopoiesis, 
and  lymphopoiesis  have  been  developed  to  better  understand  the  clinical  time  course  of  blood 
cells  after  acute  radiation  exposure.  The  models  have  been  optimized  and  validated  against 
blood  cell  kinetic  data  from  human  radiation  accidents  and  accurately  illustrate  the  average 
observed  trends  in  hematopoietic  dynamics  after  acute  radiation  exposure.  The  models  not 
only  provide  valuable  information  on  the  expected  clinical  course  after  radiation  exposure, 
but,  by  integrating  parameter  changes  for  additional  insults  such  as  thermal  injury,  they  also 
provide  a  mechanism  for  understanding  the  impact  of  combined  injuries.  By  relating  key 
clinical  metrics  obtained  from  the  models,  such  as  duration  of  neutropenia,  the  models  may 
also  assist  in  predicting  the  probability  and  time  to  mortality  from  combined  injuries,  facets 
for  which  little  to  no  data  currently  exists.  Future  work  will  further  increase  the  utility  of 
these  models  by  integrating  treatment  and  population  variability  to  support  the  application 
of  these  models  in  emergency  preparedness  scenario  analyses  and  medical  resource  planning. 
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Appendix  A  Human  Case  Study  Data 

Here  the  process  used  to  evaluate  case  study  data  for  use  in  optimization  and  validation 
of  the  mechanistic  hematopoietic  models  is  described.  Case  study  data  was  collected  from 
published  literature  on  reactor  and  industrial  accidents,  the  Chernobyl  accident,  and  early 
radiation  therapy  studies.  A  larger  number  of  case  studies  were  identified  than  are  reported 
here,  however,  many  of  them  lacked  detail  on  hematological  parameters  over  time,  had 
significant  treatment,  or  were  overtly  localized  exposures.  This  review  focuses  on  the  initial 
down-selection  of  case  studies  that  could  potentially  provide  valuable  data  for  our  work. 
Nevertheless,  many  of  the  case  studies  have  details  that  must  be  carefully  considered  before 
their  associated  data  are  used  in  our  modeling  work. 

A.l  General  Considerations  and  Criteria 

Case  studies  were  reviewed  with  consideration  of  factors  that  could  impact  hematopoietic 
response  in  humans,  which  include:  partial  body  or  non-uniform  exposures,  uncertainty 
in  dose  estimates,  differential  neutron  (or  x-ray)  energies  and  contributions  to  dose  and 
the  biological  response  to  them,  existing  individual  variability  (due  to  different  baseline  cell 
counts,  co-morbidities,  and/or  individual  responses),  and  the  application  of  treatment  and/or 
therapeutics.  Data  was  obtained  on  as  many  of  these  details  as  possible.  Some  case  studies, 
such  as  the  Marshallese  and  Chernobyl  patients,  were  exposed  to  radioactive  fallout  that 
included  beta  emitting  radionuclides  resulting  in  some  cases  to  cutaneous  injury  and  some 
level  of  internalized  contamination.  Case  studies  recommended  for  use  in  optimizations  were 
prioritized  according  to  relative  reliability  of  dose  estimates,  uniformity  of  doses,  availability 
of  time  course  data  on  the  cell  counts,  and  cases  without  significant  hematological  treatment. 

Where  possible,  doses  in  Sv  were  used  to  account  for  the  relative  biological  effectiveness 
(RBE)  of  the  neutrons  that  were  involved  in  the  criticality  accidents.  In  two  specific  cases, 
absorbed  doses  in  Gy  were  recorded  without  regard  to  neutron  RBE  due  to  uncertainties 
associated  with  the  dose  estimates  resulting  in  unrealistic  doses  given  the  observed  effects  in 
these  case  studies.  These  studies  are  discussed  below.  The  large  uncertainty  placed  on  the 
doses  estimates  for  these  case  studies  account  for  any  effect  of  neutron  RBE.  Special  care 
was  needed  to  ensure  that  doses  used  for  the  modeling  were  consistently  input  as  free-in-air 
(FIA)  or  midline  tissue  (MLT)  doses.  The  reporting  of  dose  estimates  for  the  studies  was 
highly  variable.  The  doses  reported  here  reflect  the  values  used  as  inputs  as  FIA  gamma 
equivalent  doses  for  modeling  hematopoietic  responses  across  the  different  types  of  radiation 
and  exposure  scenarios. 

Although  a  few  case  studies  that  involved  treatment  with  significant  transfusions,  cy¬ 
tokines,  and/or  stem  cell  transplantation  were  included  in  the  optimization,  data  points 
collected  after  these  treatments  were  removed  from  the  optimization  data.  Some  case  study 
data  for  treated  patients  remained  in  the  validation  data  out  of  necessity,  however,  model 
comparisons  (reported  elsewhere)  discuss  the  limitations  of  the  data  and  relevance  of  the 
comparisons. 

Our  current  modeling  work  has  the  aim  of  predicting  hematopoietic  response  resulting 
from  relatively  uniform  whole  body  exposures.  Since  significant  acute  localized  injury  re¬ 
sults  in  a  different  overall  clinical  picture  than  uniform  whole  body  exposures,  case  studies 
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involving  highly  non-uniform/partial  body  exposures  were  avoided  for  use  in  optimization. 
Such  cases  would  be  expected  to  significantly  deviate  from  our  modeling  results  as  well. 
Although  an  attempt  was  made  to  avoid  significant  partial  body  exposure  case  studies  that 
could  not  afford  reasonable  comparisons,  nearly  all  of  the  studies  involved  some  degree  of 
inhomogeneity  in  the  exposure. 

A. 2  Case  Study  Descriptions 
A. 2.1  Criticality  accidents 
Los  Alamos,  1945  and  1946 

The  criticality  accident  at  Los  Alamos  National  Laboratory  in  1945  resulted  in  two  expo¬ 
sures;  one  fatality  and  one  minor  exposure  (Hempelmann  et  al.  1952;  Hempelmann  1961; 
Hoffman  and  Hempelmann  1957;  Bond  et  al.  1965;  Mettler  et  al.  2001).  The  fatality,  referred 
to  as  Case  1,  had  sustained  severe  localized  exposure  to  the  hand,  resulting  in  significant 
complications  leading  to  death.  Case  2  had  a  relatively  minor  exposure  to  the  back.  Al¬ 
though  hematological  data  was  available  for  both  cases,  only  Case  2,  the  low  dose  exposure, 
was  used  in  the  optimization.  The  acute  partial  body  nature  of  the  exposure  in  Case  1 
prevents  its  utility  in  the  model  and  was  instead  used  to  visual  the  effects  of  non-uniform 
exposures  (5.5-20  Sv). 

Another  criticality  occurred  in  1946  in  which  eight  persons  were  exposed  (Cases  3-10). 
Hematological  data  was  available  on  seven  of  these  persons.  Cases  3  and  4  were  significant  but 
highly  non-uniform,  prohibiting  use  in  model  optimization.  Case  4  deemed  fit  for  validation, 
while  Case  3  was  used  to  visualize  the  potential  effects  of  partial  body  exposure.  Other  cases 
were  moderate  to  insignificant  exposures,  were  fairly  uniform,  and  were  not  treated.  Cases 
6-10  had  hematological  data  and  were  included  in  validation. 

It  is  important  to  note  that  the  methods  used  for  determining  platelets  for  these  studies 
differ  significantly  from  methods  applied  in  later  studies.  Therefore,  the  platelet  data  in  these 
studies  must  be  normalized  and  only  relative  platelet  changes  can  be  reliably  evaluated  in 
these  studies.  However,  several  of  the  patients  from  these  accidents  had  baseline  cell  count 
data  on  hie  since  they  were  nuclear  workers  that  were  being  monitored.  The  baseline  data 
aids  interpretation  of  the  hematological  response  after  the  exposures. 

Doses  in  Sv  as  reported  in  Mettler  et  al.  2001  were  used  for  these  case  studies  since  the 
conversion  to  Sv  accounts  for  neutron  RBE.  The  original  doses  reported  in  Hempelmann 
et  al.  1952  vary  considerably  from  the  doses  reported  later  since  neutron  doses  were  not 
accounted  for  in  the  original  work.  Also  note  that  Bond  et  al.  1965  reports  other  doses  than 
those  in  Mettler  et  al.  2001  (for  example,  Bond  reports  13.5  Gy  for  Case  3,  vs  13-21  Sv 
as  reported  in  Mettler  et  al.  2001).  However,  due  to  the  completeness  of  data  and  doses 
provided  in  Sv  found  in  Mettler  et  al.  2001,  these  are  the  data  used  for  these  case  studies. 
In  cases  where  a  dose  range  is  reported,  the  median  value  was  used  for  modeling  purposes. 

Argonne,  1952 

Four  workers  were  exposed  to  mixed  held  exposure  at  Argonne  National  Laboratory  in  1952. 
The  doses  ranged  from  1.59  to  0.11  Gy  (Bond  et  al.  1965;  Mettler  et  al.  2001).  Hematological 
data  was  obtained  on  Cases  I  and  II  (1.59  and  1.26  Gy,  respectively).  Details  on  uniformity, 
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proportion  of  neutron  contribution  to  dose,  and  treatment  were  not  available.  Since  the 
workers  were  on  a  platform  some  distance  from  the  source,  relative  uniformity  is  assumed. 
Only  Case  I  showed  any  clinical  symptoms,  however,  blood  count  changes  were  observed.  A 
median  dose  of  1.5  Gy  (1-2  Gy  range)  was  assigned  for  both  of  cases,  and  Case  I  was  used 
in  the  optimization.  There  is  some  concern  over  the  appropriate  baseline  counts  for  Case  II, 
so  this  case  was  not  included  in  the  optimization. 

Y-12,  1958 

The  criticality  accident  that  occurred  at  the  Y-12  facility  in  Oak  Ridge,  TN  in  1958  resulted 
in  eight  exposed  persons  with  significant  doses  to  five  of  the  individuals  (Mettler  et  al.  2001). 
The  doses  were  mixed  neutron  and  gamma.  Useful  hematological  data  was  obtained  on  four 
of  these  individuals:  patients  A,  B,  C,  and  E  (Mettler  et  al.  2001;  Andrews  et  al.  1961;  Bond 
et  al.  1965).  All  of  these  patients  were  in  the  room  at  the  time  of  the  criticality  at  different 
distances  from  the  source.  The  exposures  were  total  body  but  were  not  completely  uniform. 

Although  total  doses  in  Sv  were  reported  in  Mettler  et  al.  2001,  upon  review  we  found 
several  dosimetric  studies  that  reported  different  doses  for  the  case  study  (Hurst  and  Ritchie 
1959;  Hurst,  Ritchie,  and  Emerson  1959;  Mole  1984;  Kerr  and  Tankersley  2006).  In  par¬ 
ticular,  the  study  by  Mole  1984,  indicated  that  the  sodium  ratios  used  to  calculate  the 
neutron  contributions  in  the  early  dose  estimates  for  this  accident  and  the  Vinca  accident 
were  off.  Therefore,  the  revised  doses  in  Gy  as  reported  in  the  Kerr  and  Tankersley  2006 
dose  reconstruction  in  2006  were  used  as  the  MLT  dose  input  for  the  model  optimization. 
For  these  studies,  the  RBEs  were  not  accounted  for  because  the  estimates  were  unrealistic 
given  observed  health  effects  when  an  RBE  was  used  for  the  neutron  contribution.  However, 
the  uncertainty  provided  for  the  dose  estimates  is  considerable  (25%  for  both  neutron  and 
gamma  contributions). 

Although  bone  marrow  transplantation  was  considered,  patients  were  expected  to  recover 
and  in  the  end  only  received  antibiotics  as  necessary.  Due  to  the  relative  reliability  of  the 
dose  estimated,  relative  uniformity  of  exposure,  and  lack  of  treatment,  these  studies  were  all 
used  in  the  optimization  of  the  models.  The  one  exception  to  this  was  thrombopoiesis  and 
granulopoiesis  data  from  patient  E  which  was  reserved  for  validation. 

Vinca,  1958 

The  Vinca  reactor  criticality  that  occurred  in  1958  resulted  in  six  exposed  individuals  and  one 
fatality  (Mettler  et  al.  2001).  Doses  were  relatively  uniform,  however,  there  are  significant 
discrepancies  in  the  doses  reported  (Jammet  et  al.  1959;  Hurst,  Ritchie,  Sanders,  et  al.  1961; 
Pendic  1961;  Auxier  1961;  Mole  1984;  Pesic  2012).  Originally,  the  midpoint  of  the  total 
dose  range  reported  in  Jammet  et  al.  1959  was  used  for  modeling  purposes,  however,  the 
dose  ranges  were  not  realistic  given  observed  health  effects.  After  review  of  the  collective 
dosimetry  studies  for  this  accident,  the  doses  calculated  after  correcting  the  sodium  ratios 
(referenced  in  the  Y-12  accidents,  see  Mole  1984)  and  as  reported  in  Pesic  2012  were  used. 

The  patients  received  only  nominal  blood  transfusions  to  compensate  for  blood  draws. 
Bone  marrow  transplantations  were  attempted  in  five  of  the  six  cases;  therefore,  only  data 
up  to  the  point  of  transplantation  was  used.  The  majority  of  the  data  from  this  case  study 
was  obtained  from  the  original  graphs  found  in  Jammet  et  al.  1959.  However,  the  early 
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time  points  were  particularly  difficult  to  reconstruct.  Data  for  three  early  time  points  were 
reported  in  Pendic  1961  and  included  in  our  data  for  optimization. 

Los  Alamos,  1958 

A  third  criticality  at  Los  Alamos  National  Laboratory  resulted  in  an  operator  being  exposed 
to  an  extremely  high  dose  mixed  field  of  neutrons  and  gamma  radiation  (Mettler  et  al.  2001; 
Bond  et  al.  1965).  Few  data  hematological  data  points  were  available  since  the  patient  died 
at  35  hours.  Two  other  workers  were  exposed  to  1.35  and  0.54  Gy,  but  detailed  hematology 
was  not  available  on  these  cases. 

Sarov-1963 

Two  persons  (KH  and  M)  were  involved  in  this  criticality  accident,  having  received  approx¬ 
imately  5.5  and  3.7  Gy  of  neutron  and  gamma  exposure,  respectively.  The  patients  received 
blood  transfusions  and  survived  the  incident.  Data  on  patient  KH  was  obtained  from  Mettler 
et  al.  2001  and  included  for  validation. 

Mol-1965 

An  operator  standing  above  the  tank  received  a  non-uniform  dose  of  mixed  field  radiation. 
The  dose  in  the  pelvic  region  was  10  Gy  but  was  2-5  Gy  in  the  chest;  a  mean  bone  marrow 
dose  was  estimated  to  be  5  Gy.  One  heavily  irradiated  leg  (50  Gy)  was  amputated,  but  the 
patient  survived.  Details  on  hematological  was  lacking.  Due  to  the  amputation  and  highly 
localized  injury,  this  study  was  not  considered  for  optimization  or  validation.  However,  data 
was  obtained  from  Mettler  et  al.  2001  and  used  for  model  comparisons. 

Tokai-mura,  1999 

The  Japanese  criticality  resulted  in  two  highly  exposed  individuals  that  eventually  died  and 
a  third  that  was  moderately  exposed  (Hirama  et  al.  2003;  Mettler  et  al.  2001).  All  patients 
received  aggressive  care.  Data  from  subjects  A  and  B  were  included  in  the  validation,  while 
data  from  subject  C  was  compared  to  the  model  output  to  examine  potential  effects  of 
treatment. 

A. 2. 2  Radiation  therapy 

In  the  early  days  of  radiation  therapy,  a  number  of  patients  were  treated  with  single  doses  of 
whole  body  irradiation  followed  by  minimal  hematological  therapy.  The  advantage  to  data 
from  radiation  therapy  case  studies  is  that  baseline  blood  cell  counts  are  obtained  and  the 
radiation  dose  is  known  with  much  more  certainty  than  in  radiation  accidents.  The  obvious 
shortcoming  in  using  radiation  therapy  data  is  that  patients  have  co-morbidities  that  may 
impact  the  observed  hematological  effects.  From  the  literature  evaluated,  potentially  useful 
data  from  two  studies  were  identified:  one  study  on  patients  with  cancer  not  involving  the 
hematopoietic  system  and  one  study  on  patients  with  renal  dysfunction  receiving  radiation 
before  kidney  transplantation. 

Non-hematological  cancer 
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Hematological  data  was  collected  in  cancer  patients  receiving  therapeutic  exposures  to  250 
kvp  x-rays  in  three  different  dose  ranges:  2,  1.5,  and  1  Gy  (population  sizes  were  30,  12, 
and  18,  respectively).  Although  tumors  in  these  patients  originated  in  organs  systems  other 
than  the  hematopoietic  system  (kidney,  lung,  colon,  etc.),  the  disease  state  may  still  affect 
the  hematopoietic  response.  Nevertheless,  the  data  provided  good  baselines  for  each  cohort, 
relatively  precise  doses,  and  mean  and  standard  deviation  for  the  cohorts.  Unfortunately, 
the  data  only  extends  out  to  about  10  days  post-exposure. 

Renal  disorders 

Hematological  data  in  a  cohort  of  patients  receiving  radiation  therapy  before  kidney  trans¬ 
plantation  was  studied  (Dealy  and  Tubiana  1964;  Mathe  et  al.  1964;  Tubiana  et  al.  1961). 
Most  of  these  patients  received  fractionated  doses  of  250  kvp  x-rays  or  exposure  to  a  Co-60 
gamma  source.  The  fractionated  exposures  could  not  be  used  in  the  modeling  effort,  how¬ 
ever,  two  patients  in  the  cohort  received  single  4  Gy  gamma  exposures.  The  hematological 
data  was  reported  as  an  average  for  these  two  patients.  When  hemorrhages  were  observed, 
platelet  transfusions  were  performed  (Tubiana  et  ah  1961),  however,  it  was  noted  that  the 
transfusions  did  relieve  the  bleeding  but  did  not  restore  circulating  platelet  levels. 

A. 2. 3  Industrial  irradiator  accidents 
UT  CARL  Co-60,  1971 

An  individual  was  relatively  uniformly  exposed  to  Co-60  source  at  the  Variable  Dose  Rate 
Irradiation  Facility  (VDRIF;  UT  CARL,  TN  in  1971,  Goans  2012).  The  estimated  dose  was 
2.6  Gy;  details  on  treatment  were  not  available  and  are  assumed  to  be  negligible  given  the 
moderate  dose. 

Nesvizh,  Belarus  Co-60,  1991 

At  a  Nesvizh,  Belarus  commercial  irradiation  facility  in  1991,  an  individual  was  exposed  to 
a  Co-60  sterilization  source  (Baranov  et  ah  1994;  Mettler  et  ah  2001).  The  person  received  a 
relatively  uniform  whole  body  gamma  exposure  estimated  between  10-15  Gy.  The  midpoint 
of  12.5  Gy  was  used  for  modeling.  Granulocyte  and  platelet  data  was  available  for  the  case 
study,  however,  IL-3  and  granulocyte  macrophage  colony-stimulation  factor  (GM-CSF)  as 
well  as  blood  products  were  administered  from  day  3,  limiting  the  amount  of  data  available 
for  use. 

Brescia  1975 

One  individual  was  exposed  at  a  Co-60  irradiation  facility  in  Italy  with  limited  details  avail¬ 
able  (Mettler  2001).  Details  on  medical  treatment  were  not  available.  The  patient  died  after 
12  days.  This  was  used  for  validation. 

Kjeller  1982 

An  individual  was  lethally  irradiated  by  an  exposed  Co-60  source  at  a  sterilization  facility 
in  Norway  (Stavem  et  ah  1985).  The  exposure  was  non-uniform  and  there  was  considerable 
uncertainty  in  the  dose  estimate.  Antibiotic  treatments  were  delivered,  and  neutrophil, 
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platelet,  and/or  whole  blood  transfusions  were  administered  as  needed  from  7-8  days. 

Soreq  1990 

An  individual  received  a  whole  body  exposure  to  a  Co-60  source  at  an  irradiation  facility 
in  Israel  (IAEA  1993;  Mettler  et  al.  2001).  The  estimated  dose  was  between  10  and  20  Gy. 
Supportive  care  and  aggressive  therapy  were  administered.  G-CSF  was  given  on  day  3,  IL-3 
on  day  5,  platelet  transfusions  as  needed,  and  eventual  bone  marrow  transplantation.  The 
patient  died  day  36  with  a  wide  range  of  complications. 

Shanghai  1990 

A  total  of  seven  workers  were  irradiated  by  a  Co-60  source  at  an  industrial  irradiation 
facility  in  China  (Mettler  et  al.  2001).  The  exposures  were  very  non-uniform,  and  the  two 
most  heavily  irradiated  patients,  S  and  W,  received  average  whole  body  doses  of  12  and  11 
Gy,  respectively.  The  mid-section  of  the  patients  received  the  highest  doses.  Both  patients 
received  whole  blood  and  platelet  transfusions  starting  day  2.  Bone  marrow  transplantation 
was  performed  on  day  11  and  4  for  S  and  W,  respectively.  S  died  on  day  25,  and  W  died  of 
respiratory  failure  on  day  90. 

A. 2. 4  Miscellaneous 
Marshallese  1954 

The  mean  hematological  responses  observed  in  64  Marshallese  on  the  Rongelap  atoll  exposed 
to  fallout  radiation  from  the  atomic  bomb  testing  in  the  Pacific  Ocean  was  reviewed  in  Bond 
et  al.  1965.  The  average  whole-body  gamma  dose  for  this  group  was  1.68  Gy.  In  this  cohort, 
the  exposures  resulted  from  the  radioactive  fallout  and  was  considered  somewhat  uniform. 
Internal  contamination  and  skin  doses  were  not  taken  into  account  for  the  doses  reported 
here.  Some  of  the  Marshallese  did  receive  burns  to  the  skin,  and  the  doses  were  received 
over  a  period  of  time.  These  factors  limit  the  relevance  of  comparisons  to  our  acute  exposure 
hematological  models. 

Lockport  1960 

This  incident  involved  nine  persons  being  exposed  to  a  pulse  of  x-rays  from  a  powerful  radar 
generator  due  to  the  protective  shielding  having  been  removed  (Howland  et  al.  1961;  Bond 
et  al.  1965).  The  doses  received  to  the  persons  were  very  non-uniform,  ranging  from  15 
Gy  to  the  head  to  3  Gy  to  the  trunk.  Although  one  individual  was  severely  injured,  all 
nine  persons  recovered.  Data  on  two  of  the  patients  with  average  whole  body  doses  of  3 
Gy  have  been  included  in  the  validation  for  comparison;  however,  Bond  et  al.  1965  notes 
the  unusually  slow  recovery  as  compared  to  uniform  whole  body  doses.  The  RBE  of  the 
x-rays,  the  non- uniformity,  and  the  uncertainty  in  the  doses  may  contribute  to  the  unusual 
observations  in  this  study. 

Chernobyl  1986 

The  average  hematological  response  of  eleven  patients  Chernobyl  nuclear  power  plant  ac¬ 
cident  was  summarized  and  reviewed  by  Fliedner,  Grade,  et  al.  2007.  This  cohort  was 
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referenced  as  having  level  3  hematopoietic  syndrome  with  an  estimated  whole  body  dose  of 
4.85  Gy  gamma  dose.  It  was  noted  that  these  patients  experienced  hematopoietic  recov¬ 
ery  without  cytokine  therapy.  This  cohort  includes  beta  irradiation,  skin  contamination, 
and  internal  contamination.  Exposures  were  obtained  over  an  extended  period  of  time  that 
spanned  up  to  several  weeks. 

UNSCEAR  1988:  Case  Studies  from  Chernobyl 

The  overall  experience  from  the  Chernobyl  accident  was  reviewed  in  the  Guskova  et  al.  1988 
report.  In  this  review,  the  hematological  responses  of  four  moderately  exposed  patients  are 
provided.  It  was  noted  that  treatment  for  these  patients  was  conservative,  with  no  antibiotics 
or  transfusions  unless  symptoms  warranted.  However,  no  specific  mention  of  administration 
of  either  of  these  was  provided.  One  patient  in  particular  was  noted  to  have  significant  skin 
burns  (case  39).  These  patients  received  no  specialized  therapy.  As  indicated  previously,  the 
Chernobyl  cohort  includes  beta  irradiation,  skin  contamination,  internal  contamination,  and 
protracted  exposures. 

Russian  Examples 

Exposures  occurring  in  the  former  Soviet  Union  described  originally  by  Guskova  and  Baranov 
in  Russian,  were  cited  for  illustrative  examples  for  a  radiobiological  review  of  radiation 
sickness  in  man  (Yarmonenko  1988;  Baranov  1981).  Therefore,  details  of  these  exposures 
are  vague  but  note  that  they  are  both  non-uniform,  neutron  and  gamma  exposures.  Details 
on  treatment  are  lacking. 

A. 3  Summary  and  Conclusions 

A  sufficient  amount  of  data  for  optimization  and  validation  of  the  hematopoietic  model  under 
development  was  obtained  from  peer  review  literature  on  radiation  accident  case  studies  and 
therapy  patients.  However,  much  care  was  needed  in  consistently  assigning  radiation  doses 
and  controlling  for  therapeutic  factors  in  observed  hematopoietic  responses.  Although  much 
of  the  data  was  removed  due  to  treatment  and/or  significant  partial  body  effects,  these 
additional  data  or  case  studies  can  be  used  in  the  future  to  better  understand  the  impact 
of  partial  body  exposures  and  bone  marrow  sparring  as  well  as  the  impact  of  treatment  on 
hematopoietic  recovery. 
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Table  A.l:  (1  of  2)  Case  studies  used  in  optimization 


Case  Study-Patient 

References 

Estimated  Dose;  Range 
(Type) 

Exposure  Details 

Treatment  Details 

Los  Alamos  1945-2obc 

Hempelmann  et  al.  1952 
Hempelmann  1961 

Hoffman  and  Hempelmann  1957 
Mettler  et  al.  2001 

0.12  Gy  (n,  7) 

Insignificant  exposure  to 
back 

None 

Argonne  1952-Iabc 

Bond  et  al.  1965 

Mettler  et  al.  2001 

1 .59;  1-2  Gy  (n,  7) 

Assumed  uniform 

Some  clinical  effects,  assume  no 
significant  treatment 

Y-12  1958-Aabc 

Bond  et  al.  1965 

Mettler  et  al.  2001 

Andrews  et  al.  1961 

Kerr  and  Tankersley  2006 

4.02  Gy±25%  (n,  7) 

Criticality,  relatively  uni¬ 
form  exposure 

Antibiotics 

Y-12  1958-Bobc 

II 

2.97  Gy±25%  (n,  7) 

II 

II 

Y-12  1958-Cabc 

II 

3.73  Gy±25%  (n,  7) 

II 

None 

Y-12  1958-EC 

II 

2.59  Gy±25%  (n,  7) 

II 

II 

Vinca  1958-Vabc 

Jammet  et  al.  1959 

Jammet  1961 

Pendic  1961 

Pesic  2012 

6.11  Gy±15%  (n,  7) 

Criticality,  relatively  uni¬ 
form  exposure 

Bone  marrow  transfusion  wk  5 
small  blood  transfusion 
platelets  during  hemorrhage 
supportive  care 

Vinca  1958-Dabc 

II 

5.4  Gy±15%  (n,  7) 

II 

II 

Vinca  1958-HQbc 

II 

4.37  Gy±15%  (n,  7) 

II 

II 

Vinca  1958-Mabc 

II 

5.74  Gy±15%  (n,  7) 

II 

Bone  marrow  transfusion  wk  5 
small  blood  transfusion 
supportive  care 

Vinca  1958-Gabc 

II 

5.35  Gy±15%  (n,  7) 

II 

II 

Vinca  1958-Babc 

II 

3.51  Gy±15%  (n,  7) 

II 

Small  blood  transfusion 
supportive  care 

Cancer  patients  1958- 
mean  of  30abc 

Miller  et  al.  1958 

2  Gy  (250  kvp  x-ray) 

Uniform,  whole-body 

None 

Cancer  patients  1958- 
mean  of  18abc 

II 

1  Gy  (250  kvp  x-rays) 

II 

II 

Cancer  patients  1958- 
mean  of  12abc 

II 

1.5;1.25-1.75  Gy 
(250  kvp  x-rays) 

II 

II 

a  Case  study  used  in  thrombopoiesis  optimization 
b  Case  study  used  in  granulopoiesis  optimization 
c  Case  study  used  in  lymphopoiesis  optimization 


Table  A.l:  (2  of  2)  Case  studies  used  in  optimization 


Case  Study-Patient 

References 

Estimated  Dose;  Range 
(Type) 

Exposure  Details 

Treatment  Details 

UT  CARL  1971“bc 

Goans  2012 

2.6  Gy  (7) 

Relatively  uniform 

No  treatment  assumed 

Nesvizh  1991“ 

Baranov  et  al.  1994 

12.5;10-15  Gy  (7) 

Relatively  uniform 

Cytokines  and  blood  transfusions 

Lockport  1960-AC 

Bond  et  al.  1965 

Howland  et  al.  1961 

3  Gy 
-x-rays 

Very  non-uniform 

Conservative  treatment  noted 
Antibiotics  and  transfusions  only  if 
necessary 

Lockport  1960-BC 

II 

11 

II 

II 

Yarmonenko  1988- 

r^abc 

II 

9.8  Sv 
-n,  7 

Non-uniform 

No  details 

a  Case  study  used  in  thrombopoiesis  optimization 
b  Case  study  used  in  granulopoiesis  optimization 
c  Case  study  used  in  lymphopoiesis  optimization 


Table  A. 2:  (1  of  3)  Case  Studies  for  potential  use  in  validation  and  future  work 


Case  Study-Patient 

References 

Estimated  Dose; 
Range  (Type) 

Exposure  Details 

Treatment  Details 

Notes  on  applicability 

Most  reliable  case  studies 

Los  Alamos  1946-4abc 

Los  Alamos  1946-6abc 

Los  Alamos  1946-7abc 

Los  Alamos  1946-8abc 

Los  Alamos  1946-9abc 

Los  Alamos  1946-10abc 

Hempelmann  et  al.  1952 
Hempelmann  1961 

Hoffman  and  Hempelmann 
1957 

Mettler  et  al.  2001 

II 

II 

II 

II 

II 

3.6  Sv  (n,  7) 

1.6  Sv  (n,  7) 

1.1  Sv  (n,  7) 

0.65  Sv  (n,  7) 
0.47  Sv  (n,  7) 
0.37  Sv  (n,  7) 

Non-uniform  criticality, 
some  localized  exposure 

Relative  uniform  criti¬ 
cality 

II 

II 

II 

II 

Antibiotics,  fluids,  blood 
transfusions 

No  effects  observed,  no 
treatment 

II 

II 

II 

II 

Limitations: 

some  treatment 

Argonne  1952-IIC 

Bond  et  al.  1965 

Mettler  et  al.  2001 

1.26;  1-2  Gy  (n, 
7) 

Assumed  uniform 

No  clinical  effects,  no  treat¬ 
ment  assumed 

Y-12  1958-Eab 

Bond  et  al.  1965 

Mettler  et  al.  2001 

Andrews  et  al.  1961 

Kerr  and  Tankersley  2006 

2.59  Gy±25% 

(n,  7) 

Criticality,  relatively 
uniform  exposure 

Antibiotics 

Lockport  1960-Aab 

Bond  et  al.  1965 

Howland  et  al.  1961 

3  Gy  (x-rays) 

Very  non-uniform 

Conservative  treatment 

noted,  antibiotics  and 

transfusions  only  if  neces¬ 
sary 

Limitations: 

non-uniform 

Lockport  1960-Bab 

II 

II 

II 

II 

II 

Kidney  transplants  1964 
-mean  of  2d 

Mathe  et  al.  1964 

Tubiana  et  al.  1961 

Dealy  and  Tubiana  1964 

4  Gy  (7) 

Uniform,  whole-body 

Platelet  transfusions  noted 
amelioration  of  hemorrhage 
but  not  platelet  levels 

Limitations: 
repeated  transfusion 

Kjeller  1982abc 

Stavern  et  al.  1985 

Mettler  2001 

22.5  Gy  (7) 

Dosimetry  uncertain, 
non-uniform 

Antibiotics,  transfusions  on 
days  8,  9,  12;  died  day  13 

Limitations:  non- 

uniform  and  treatment 

a  Case  study  used  in  thrombopoiesis  validation  c  Case  study  used  in  lymphopoiesis  validation 

b  Case  study  used  in  granulopoiesis  validation  d  Case  study  used  to  visualize  treatment  effects 


Table  A. 2:  (2  of  3)  Case  studies  for  potential  use  in  validation  and  future  work 


Case  Study-Patient 

References 

Estimated  Dose; 
Range  (Type) 

Exposure  Details 

Treatment  Details 

Notes  on  applicability 

Case  studies  with  mild  limitations 

Marshallese 
mean  of  64abc 

1954- 

Bond  et  al.  1965 

1.68  Gy  (7,  P) 

Uniform  fallout,  in¬ 
cluded  skin  contamina¬ 
tion 

No  details,  none  assumed 

Limitations:  cutaneous  doses, 
internal  contamination  and 
protracted  doses 

Sarov  1963abc 

Mettler  et  al.  2001 

5.5  Gy  (n,  7) 

Non-uniform  criticality 

Antibiotics  and  blood 
transfusions  every  3-5  days 

Limitations:  non-uniform  and 

treatment 

Brescia  1975ab 

Mettler  2001 

12  Gy  (7) 

No  details 

No  details,  died  day  12 

Limitations:  details  on  expo¬ 
sure  and  treatment  lacking 

Chernobyl 
mean  of  11 

1986- 

Fliedner,  Grafile,  et  al. 
2007 

4.85;2.6-7.1  Gy 

(7)  b) 

Variable  exposure,  may 
include  beta  contami¬ 
nation 

No  cytokine  therapy, 

hematopoietic  recovery 

Limitations:  large  dose  range, 
cutaneous  doses,  internal 
contamination  and  protracted 
doses 

Chernobyl 
Case  21 

1986- 

Guskova  et  al.  1988 

3.9;3.6-4.2  Gy 
(7>  /3) 

Variable  exposure,  beta 
and  internal  contami¬ 
nation 

Supportive  care  only 

Limitations:  skin  burns,  in¬ 
ternal  contamination  and  pro¬ 
tracted  doses 

Chernobyl 
Case  39 

1986- 

II 

2.85;2.4-3.3  Gy 
(7:  b) 

II 

II 

Limitations:  cutaneous  doses, 
internal  contamination  and 
protracted  doses 

Chernobyl 
Case  48 

1986- 

II 

1.25;1. 1-1.4  Gy 

(7,  b) 

II 

II 

II 

Chernobyl 
Case  97 

1986- 

II 

0.6;0.3-0.9  Gy 

(7)  /3) 

II 

II 

II 

Shanghai  1990-Sabc  Mettler  2001 

12  Gy  (7) 

Very  non-uniform 

Antibiotics,  whole  blood 
and  platelet  transfusions, 
bone  marrow  transplant 

Limitations:  details  on  non¬ 
uniformity  not  available  and 
significant  treatment 

Shanghai  1990- Wd  " 

11  Gy  (7) 

II 

II 

Limitations:  some  non¬ 

uniformity  and  significant 
treatment 

a  Case  study  used  in  thrombopoiesis  validation 
b  Case  study  used  in  granulopoiesis  validation 


c  Case  study  used  in  lymphopoiesis  validation 
d  Case  study  used  to  visualize  treatment  effects 


Table  A. 2:  (3  of  3)  Case  studies  for  potential  use  in  validation  and  future  work 


Case  Study-Patient 

References 

Estimated  Dose; 
Range  (Type) 

Exposure  Details 

Treatment  Details 

Notes  on  applicability 

Tokai-mura  1999-Aabc 

Hirama  et  al.  2003 

>20  GyEq  (n, 

Non-uniform  criticality 

Antibiotics,  G-CSF  day  1, 

Limitations:  non-uniform 

Mettler  et  al.  2001 

7) 

eventual  stem  cell  trans¬ 
plant,  died  after  82  days 

and  some  treatment 

Tokai-mura  1999-Babc 

II 

7.4;6-10  GyEq 

Non-uniform  criticality 

Antibiotics,  G-CSF  day  1, 

II 

(n>  7) 

eventual  stem  cell  trans¬ 
plant,  died  after  210  days 

Tokai-mura  1999-Cd 

II 

2.3;  1-4.5  GyEq 

Relatively  uniform  criti- 

Antibiotics,  G-CSF  day  2, 

Limitations:  treatment 

(n>  7) 

cality 

recovered 

Case  studies  with  dramatic  limitations 

Los  Alamos  1945-1® 

Hempelmann  et  al.  1952 

12.75;5.5-20  Sv 

Non-uniform  criticality, 

Antibiotics,  died  day  25 

Not  useful:  Acute  localized 

Hempelmann  1961 

(n»  7) 

acute  localized  exposure 

exposure;  useful  in  showing 

Hoffman  and  Hempel¬ 
mann  1957 

Mettler  et  al.  2001 

partial  body  effects 

Los  Alamos  1946-3® 

II 

17;  13-20  Sv  (-n, 

Non-uniform  criticality, 

Antibiotics,  fluids,  blood 

Not  useful:  Limited  data 

7) 

acute  localized  exposure 

transfusions,  died  day  9 

and  acute  localized  expo- 

sure 

Mol  1965ab® 

Mettler  et  al.  2001 

5.5  Gy  (n,  7) 

Non-uniform  criticality 

Antibiotics  and  arnputa- 

Not  useful:  Acute  localized 

tion 

exposure  resulting  in  ampu¬ 
tation 

Baranov  1981 

Baranov  1981 

2.3;2-2.6  Gy  (7) 

No  details 

No  details,  assume  no 

Not  useful:  no  details 

treatment 

available  on  exposure,  data 
inconsistent  with  known 
hematological  response 

Yarmonenko  1988-K® 

Yarmonenko  1988 

5.8;2.8-10  Sv  (n, 

Non-uniform:  10  and  2.8 

No  stem  cell  transplant, 

Not  useful:  Highly  non- 

7) 

Sv  to  each  side 

details  limited,  sponta- 

uniform;  useful  in  showing 

neous  recovery 

partial  body  effects 

Soreq  1990d 

IAEA  1993 

15;  10-20  Gy  (7) 

Whole-body 

Antibiotics,  G-CSF  day 

Not  useful:  Highly  treated; 

Mettler  2001 

1,  IL-3  day  5,  platelet 

useful  in  showing  effect  of 

transfusions,  bone  marrow 
transplant,  died  day  36 

treatment 

a  Case  study  used  in  thrombopoiesis  validation  c  Case  study  used  in  lymphopoiesis  validation  e  Case  study  used  to  visualize  partial  body  effects 

b  Case  study  used  in  granulopoiesis  validation  d  Case  study  used  to  visualize  treatment  effects 


Appendix  B  Supplemental  Model  Comparisons 

B.l  Optimization  Data  Comparisons 

This  section  provides  all  the  comparisons  between  the  optimization  data  and  the  correspond¬ 
ing  model  simulation  not  included  in  the  the  main  body  of  the  report.  Figure  B.l,  Figure 
B.2,  and  Figure  B.3  show  these  comparisons  for  thrombopoiesis,  granulopoiesis,  and  lym¬ 
phopoiesis,  respectively.  Each  dataset  presented  in  these  plots  conies  from  the  specified  case 
study /patient.  Details  on  each  optimization  case  study  are  provided  in  Appendix  A. 
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Figure  B.l:  (1  of  4)  Thrombopoiesis  model  compared  to  optimization  data. 
Model  output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose 
range  is  provided,  a  shaded  region. 
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(e)  Dose:  3.51  Gy  ±  15% 


(d)  Dose:  2.6  Gy 


(f)  Dose:  3.73  Gy  ±  25% 


Figure  B.l:  (2  of  4)  Thrombopoiesis  model  compared  to  optimization  data. 
Model  output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose 
range  is  provided,  a  shaded  region. 


Platelet  concentration  (dimensionless)  Platelet  concentration  (dimensionless) 


(g)  Dose:  4.37  Gy  ±  15% 


(h)  Dose:  5.35  Gy  ±  15% 


(i)  Dose:  5.4  Gy  ±  15% 


(j)  Dose:  5.74  Gy  ±  15% 


Figure  B.l:  (3  of  4)  Thrombopoiesis  model  compared  to  optimization  data. 
Model  output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose 
range  is  provided,  a  shaded  region. 


Time  following  radiation  (days) 


(k)  Dose:  6.11  Gy  ±  15% 


(1)  Dose:  12.5  Gy 


Figure  B.l:  (4  of  4)  Thrombopoiesis  model  compared  to  optimization  data. 
Model  output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose 
range  is  provided,  a  shaded  region. 
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Figure  B.2:  (1  of  3)  Granulopoiesis  model  compared  to  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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CO 


(c)  Dose:  2.0  Gy 


CO 


Time  following  radiation  (days) 

(d)  Dose:  2.6  Gy 


(e)  Dose:  3.73  Gy  ±  25% 


CO 


(f)  Dose:  4.37  Gy  ±  15% 


Figure  B.2:  (2  of  3)  Granulopoiesis  model  compared  to  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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(g)  Dose:  5.4  Gy  ±  15% 


(h)  Dose:  5.35  Gy  ±  15% 


(/) 


U) 


(i)  Dose:  5.74  Gy  ±  15%  (j)  Dose:  6.11  Gy  ±  15% 

Figure  B.2:  (3  of  3)  Granulopoiesis  model  compared  to  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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Figure  B.3:  (1  of  3)  Lymphopoiesis  model  compared  to  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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(g)  Dose:  3.73  Gy  ±  25%  (h)  Dose:  4.02  Gy  ±  25% 

Figure  B.3:  (2  of  3)  Lymphopoiesis  model  compared  to  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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C/)  w 


(i)  Dose:  5.35  Gy  ±  15% 


(j)  Dose:  5.4  Gy  ±  15% 


(k)  Dose:  5.74  Gy  ±  15% 


(1)  Dose:  6.11  Gy  ±  15% 


Figure  B.3:  (3  of  3)  Lymphopoiesis  model  compared  to  optimization  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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B.2  Other  Model  Validations  and  Comparisons 

This  section  presents  comparisons  of  the  hematopoietic  models  to  blood  cell  kinetic  data 
obtained  on  radiation  accident  victims  not  shown  in  the  main  body  of  the  report  and  not 
included  in  optimization.  Some  of  these  comparisons  were  used  for  validation,  while  others 
are  shown  to  explore  the  effects  of  various  factors  on  hematopoiesis.  Again,  details  on  each 
case  study  are  given  in  Appendix  A. 

Model  validation  comparisons  are  shown  in  Figure  B.4,  Figure  B.5,  and  Figure  B.6  for 
thrombopoiesis,  granulopoiesis,  and  lymphopoiesis,  respectively.  Studies  presented  in  this 
figure  were  considered  adequate  for  validation  purposes.  An  exception  to  this  is  the  Mol  1965 
case  study  for  thrombopoiesis  and  granulopoiesis  which  was  a  highly  non-uniform  exposure. 
The  Mol  exposure  is  included  in  the  validation  figures  to  provide  a  comparison  with  the 
Sarov  1963  case  study,  which  is  intended  for  validation,  who  received  the  same  estimated 
dose  of  5.5  Gy.  Below  we  describe  and  present  other  comparisons  not  used  for  validation. 
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Figure  B.4:  (1  of  3)  Thrombopoiesis  model  compared  to  validation  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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(c)  Dose:  0.65  Gy 


(e)  Dose:  1.6  Gy 


(d)  Dose:  1.1  Gy 


(f)  Dose:  3.0  Gy 


Figure  B.4:  (2  of  3)  Thrombopoiesis  model  compared  to  validation  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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Platelet  concentration  (dimensionless)  Platelet  concentration  (dimensionless) 


(g)  Dose:  5.5  Gy 


(i)  Dose:  20.0  Gy 


(h)  Dose:  7.4  Gy  (6-10  Gy) 


(j)  Dose:  22.5  Gy 


Figure  B.4:  (3  of  3)  Thrombopoiesis  model  compared  to  validation  data.  Model 
output  at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is 
provided,  a  shaded  region. 
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Time  following  radiation  (days) 

(b)  Dose:  0.47  Gy 


C/5 


(d)  Dose:  1.6  Gy 


Figure  B.5:  (1  of  3)  Granulopoiesis  model  compared  to  validation  data.  Output 
at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is  provided, 
a  shaded  region.  If  a  generic  baseline  value  was  used,  output  initialized  at  the 
upper  and  lower  end  of  the  healthy  granulocyte  range  is  shown  by  dotted  lines. 
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Granulocyte  concentration  (dimensionless)  Granulocyte  concentration  (dimensionles: 


(e)  Dose:  1.68  Gy 


(g)  Dose:  5.5  Gy 


C/5 


(f)  Dose:  3.0  Gy 


C/5 


Time  following  radiation  (days) 


(h)  Dose:  7.4  Gy 


Figure  B.5:  (2  of  3)  Granulopoiesis  model  compared  to  validation  data.  Output 
at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is  provided, 
a  shaded  region.  If  a  generic  baseline  value  was  used,  output  initialized  at  the 
upper  and  lower  end  of  the  healthy  granulocyte  range  is  shown  by  dotted  lines. 
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Time  following  radiation  (days) 


Time  following  radiation  (days) 


(i)  Dose:  20  Gy  (j)  Dose:  22.5  Gy 

Figure  B.5:  (3  of  3)  Granulopoiesis  model  compared  to  validation  data.  Output 
at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is  provided, 
a  shaded  region.  If  a  generic  baseline  value  was  used,  output  initialized  at  the 
upper  and  lower  end  of  the  healthy  granulocyte  range  is  shown  by  dotted  lines. 


Time  following  radiation  (days)  Time  following  radiation  (days) 

(a)  Dose:  0.37  Gy  (b)  Dose:  0.47  Gy 

Figure  B.6:  (1  of  3)  Lymphopoiesis  model  compared  to  validation  data.  Output 
at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is  provided, 
a  shaded  region.  Output  initialized  at  the  upper  and  lower  end  of  the  healthy 
lymphocyte  range  is  shown  by  dotted  lines. 
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(c)  Dose:  1.1  Gy 


(e)  Dose:  1.6  Gy 


(d)  Dose:  1.5  Gy  (1-2  Gy) 


(f)  Dose:  1.68  Gy 


Figure  B.6:  (2  of  3)  Lymphopoiesis  model  compared  to  validation  data.  Output 
at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is  provided, 
a  shaded  region.  Output  initialized  at  the  upper  and  lower  end  of  the  healthy 
lymphocyte  range  is  shown  by  dotted  lines. 
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CO 


co 


Time  following  radiation  (days) 


(g)  Dose:  7.4  Gy 


(h)  Dose:  12.0  Gy 


CO 


Time  following  radiation  (days) 


(i)  Dose:  20  Gy 

Figure  B.6:  (3  of  3)  Lymphopoiesis  model  compared  to  validation  data.  Output 
at  the  specified  dose  is  delineated  by  a  black  line  and,  if  a  dose  range  is  provided, 
a  shaded  region.  Output  initialized  at  the  upper  and  lower  end  of  the  healthy 
lymphocyte  range  is  shown  by  dotted  lines. 
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Skin  and  internal  contamination  and  protracted  exposure 

The  Chernobyl  accident  resulted  in  exposures  that  had  both  skin  exposure,  internal  contam¬ 
ination  and  were  protracted.  Case  studies  from  the  Chernobyl  accident  demonstrate  a  much 
more  severe  hematopoietic  response  than  what  our  models  predict.  Figure  B.7  shows  the 
mean  platelet  and  granulocyte  response  of  11  subjects  from  the  Chernobyl  incident  (Flied- 
ner,  GraBle,  et  al.  2007)  compared  with  the  model  output  using  the  mean  dose  value  of  4.85 
Gy.  Both  the  granulocytes  and  platelets  show  a  faster  drop  and  deeper  nadir  than  the  model 
predicts. 

Figure  B.8,  Figure  B.9,  and  Figure  B.10  show  individual  responses  of  Chernobyl  victims 
for  the  platelet,  granulocyte,  and  lymphocyte  response,  respectively.  Again,  for  both  the 
platelets  and  granulocytes  there  is  a  stronger  response  in  the  data  than  in  the  model.  How¬ 
ever,  for  the  lymphocytes  the  data  is  closer  to  what  the  model  predicts.  Radiation  dose  from 
skin  contamination  can  result  in  cutaneous  injury  leading  to  systemic  responses  and  may  im¬ 
pact  the  observed  hematopoietic  response  in  the  Chernobyl  workers.  Internal  contamination 
could  also  contribute  to  blood  cell  depletion,  and  the  protracted  nature  of  these  exposures 
further  complicate  the  hematopoietic  response. 

7/7 


Figure  B.7:  Mean  hematopoietic  response  of  11  subjects  from  Chernobyl  1986 
who  received  a  mean  dose  of  4.85  Gy  (2. 6-7.1  Gy).  Model  output  at  the  dose 
(range)  is  delineated  by  a  black  line  (shaded  region).  If  a  generic  baseline  value 
was  used,  output  initialized  at  the  upper  and  lower  end  of  the  healthy  blood  cell 
range  is  shown  by  dotted  lines. 
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Platelet  concentration  (dimensionless)  Platelet  concentration  (dimensionless) 


(a)  Dose:  0.6  Gy  (0.3-0. 9  Gy) 


(b)  Dose:  1.25  Gy  (1. 1-1.4  Gy) 


(c)  Dose:  2.85  Gy  (2.4-3. 3  Gy)  (d)  Dose:  3.9  Gy  (3. 6-4. 2  Gy) 

Figure  B.8:  Chernobyl  case  studies:  Platelet  data.  Model  output  at  the  specified 
dose  (range)  is  delineated  by  a  black  line  (shaded  region). 
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(a)  Dose:  0.6  Gy  (0.3-0. 9  Gy) 


C/5 


(b)  Dose:  1.25  Gy  (1. 1-1.4  Gy) 


C/5  C/5 


(c)  Dose:  2.85  Gy  (2.4-3. 3  Gy) 


(d)  Dose:  3.9  Gy  (3.6-4.2  Gy) 


Figure  B.9:  Chernobyl  case  studies:  Granulocyte  data.  Model  output  at  the 
specified  dose  (range)  is  delineated  by  a  black  line  (shaded  region).  Output 
initialized  at  the  upper  and  lower  end  of  the  healthy  granulocyte  range  is  shown 
by  dotted  lines. 
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(a)  Dose:  0.6  Gy  (0.3-0. 9  Gy) 


CO 


(b)  Dose:  1.25  Gy  (1. 1-1.4  Gy) 


(c)  Dose:  2.85  Gy  (2.4-3. 3  Gy) 


(d)  Dose:  3.9  Gy  (3.6-4.2  Gy) 


Figure  B.10:  Chernobyl  case  studies:  Lymphocyte  data.  Model  output  at  the 
specified  dose  (range)  is  delineated  by  a  black  line  (shaded  region).  Output 
initialized  at  the  upper  and  lower  end  of  the  healthy  lymphocyte  range  is  shown 
by  dotted  lines. 
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Time  following  radiation  (days) 


Time  following  radiation  (days) 


Time  following  radiation  (days) 

Figure  B.ll:  Example  of  variability  between  the  hematopoietic  response  of  two 
radiation  victims  receiving  2.3  Gy  (dose  range  shown  for  Tokai— mura  C).  Output 
at  2.3  Gy  and  between  1.0-4. 5  Gy  is  delineated  by  a  black  line  and  shaded  region. 
If  a  generic  baseline  value  was  used,  output  initialized  at  the  upper  and  lower 
end  of  the  healthy  blood  cell  range  is  shown  by  dotted  lines. 


Dose  and  individual  discrepancies 

Here  we  compare  two  case  studies,  Baranov  1981  and  Takai-mura  1999  C,  that  both  received 
an  estimated  2.3  Gy  radiation  dose.  Figure  B.ll  shows  the  blood  cell  kinetics  response  for 
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these  two  subjects  as  well  as  the  predicted  model  output  for  platelets,  granulocytes,  and 
lymphocytes.  For  all  three  blood  cell  types,  the  Baranov  1981  subject  had  a  more  severe 
hematopoietic  response.  This  suggests  that  the  differences  in  blood  cell  kinetics  between 
these  two  studies  is  an  issue  with  the  calculated  dose,  but  the  discrepancies  may  also  be  due 
to  treatment  effects  and  individual  variability.  Few  details  are  available  on  the  exposure  in 
the  Baranov  case  study  and  the  Tokai-mura  1999-C  study  was  highly  treated.  The  treatment 
may  partially  explain  why  the  blood  cell  kinetic  response  of  the  Tokai-mura  subject  C  is  less 
severe  than  the  Baranov  case  study. 

Treatment 

Figures  B.12  and  B.13  show  blood  cell  kinetics  from  the  kidney  transplants  study  and  the 
Shanghai  1990  W  case  study,  respectively.  Both  of  these  subjects  were  highly  treated.  For 
the  thrombopoiesis  and  granulopoiesis  systems,  the  kidney  transplant  study  has  a  faster 
hematopoietic  recovery  than  the  model  predicts.  The  kidney  transplant  subjects  received 
platelet  transfusions  which  may  have  attributed  to  this  faster  recovery.  The  granulocyte 
nadir  in  the  kidney  transplant  data  is  actually  lower  and  earlier  than  the  model  predicts 
which  is  not  what  we  would  expect  with  treatment.  On  the  other  hand,  the  disease  state  of 
the  individuals  may  impact  the  hematopoietic  response. 

The  Shanghai  W  patient  received  whole  blood,  platelet,  and  bone  marrow  transfusions, 
however  the  effect  of  this  treatment  is  not  clear  due  to  the  limited  blood  cell  data  available. 
For  thrombopoiesis,  it  appears  as  through  the  system  takes  longer  to  decline  than  the  model 


Time  following  radiation  (days) 


Figure  B.12:  Effects  of  treatment:  4.0  Gy  case  study.  Model  output  at  the 
specified  dose  is  delineated  by  a  black  line. 


109 


C/5 


Time  following  radiation  (days) 


C/5 


Time  following  radiation  (days) 

Figure  B.13:  Effects  of  treatment:  11.0  Gy  case  study.  Model  output  at  11.0 
Gy  is  delineated  by  a  black  line.  If  a  generic  baseline  value  was  used,  output 
initialized  at  the  upper  and  lower  end  of  the  healthy  blood  cell  range  range  is 
shown  by  dotted  lines. 


predicts.  However,  granulocyte  and  lymphocyte  concentrations  do  not  drastically  differ 
from  the  model  predictions.  Again,  these  deviations  may  be  due  to  both  dose  and  individual 
variabilities. 
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C/5 


Time  following  radiation  (days) 

Figure  B.14:  Effects  of  non-uniform  exposure:  5.8  Gy  (2.8-10  Gy)  case  study. 
Model  output  at  the  specified  dose  and  dose  range  is  delineated  by  a  black  line 
and  shaded  region.  If  a  generic  baseline  value  was  used,  output  initialized  at  the 
upper  and  lower  end  of  the  healthy  blood  cell  range  is  shown  by  dotted  lines. 


Non-uniform  exposures 

Figure  B.14  and  B.15  show  blood  cell  kinetics  from  non-uniform  case  studies,  including 
subject  K  from  Yarmonenko  1988  and  Los  Alamos  1946  subject  3.  Subject  K  received  10 
Gy  to  one  side  of  the  body  and  2.8  Gy  to  the  other  side,  and  model  output  is  shown  across 
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Time  following  radiation  (days) 


Time  following  radiation  (days) 


C/5 


Time  following  radiation  (days) 

Figure  B.15:  Effects  of  non-uniform  exposure:  17  Gy  (13-20  Gy)  case  study. 
Model  output  at  the  specified  dose  and  dose  range  is  delineated  by  a  black  line 
and  shaded  region.  If  a  generic  baseline  value  was  used,  output  initialized  at  the 
upper  and  lower  end  of  the  healthy  blood  cell  range  is  shown  by  dotted  lines. 


this  entire  range;  yet  the  data  still  does  not  match  what  the  model  predicts.  This  suggests 
other  factors  affect  the  response  in  this  case.  The  Los  Alamos  subject  3  also  suffered  a  very 
non-uniform  exposure;  however,  due  to  the  limited  blood  cell  data  available  it  is  difficult  to 
assess  how  this  affected  the  blood  cell  kinetics.  The  platelets  appear  to  plateau  more  prior 
to  declining,  while  the  lymphocytes  decline  sooner  and  further  than  the  model  predicts. 
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The  implication  is  that  non-uniform  exposures  may  result  in  a  more  complex  hematopoietic 
response  than  can  be  accounted  for  on  the  basis  of  radiation  dose  alone.  This  is  due  to  more 
systemic  responses  from  the  acutely  irradiated  areas  of  the  body  complicating  the  overall 
hematopoietic  response. 
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Appendix  C  Markov  Chain  Monte  Carlo  and  Identifi- 

ability  Analysis 

This  section  presents  a  more  in-depth  analysis  of  the  Markov  Chain  Monte  Carlo  (MCMC) 
and  identifiability  analysis  performed  on  the  three  optimized  hematopoietic  models.  The  goal 
of  the  MCMC  analysis  is  to  provide  insight  on  parameter  distributions  and  obtain  parameter 
confidence  intervals.  This  is  done  by  testing  different  parameter  combinations  to  see  which 
combinations  result  in  a  satisfactory  model  output.  An  identifiability  analysis  provides 
insight  into  the  quality  of  the  optimization  and  whether  the  model  is  overparameterized. 
Although  satisfactory  visual  results  can  be  obtained  with  an  overparameterized  model,  the 
lack  of  identifiability  implies  a  lack  of  biological  meaning  behind  the  obtained  parameter 
values.  For  a  more  complete  description  of  this  procedure  see  Section  2. 

Following  a  MCMC  analysis,  the  obtained  parameter  distributions  are  shown  in  Figures 
C.l,  C.2,  and  C.3  for  thrombopoiesis,  granulopoiesis,  and  lymphopoiesis,  respectively.  Ideally 
these  parameter  distributions  should  be  unimodal  with  a  clear  peak  value.  If  the  peak  occurs 
at  the  upper  or  lower  end  of  the  bounding  values,  it  implies  that,  if  biologically  reasonable, 
the  boundary  values  should  be  changed. 

Figures  C.l,  C.2,  and  C.3  also  provide  quantitative  information  on  the  correlation  be¬ 
tween  parameters.  Parameters  are  highly  correlated  if  changing  one  parameter  can  make  up 
for  another  changed  parameter  to  obtain  the  same  quality  of  optimization.  High  correlation 
exists  when  the  correlation  approaches  -1  (negatively  correlated)  or  1  (positively  correlated). 
When  two  parameters  are  highly  correlated,  obtaining  a  biologically  meaningful  optimiza¬ 
tion  is  nearly  impossible.  In  the  thrombopoiesis  model,  the  optimized  parameters  7  and  Hf) 
have  the  largest  correlation  which  equals  0.79.  In  the  granulopoiesis  model  optimization, 
the  strongest  correlation  exists  between  0  and  l  with  a  value  of  -0.82.  In  the  lymphopoiesis 
model  optimization,  the  strongest  correlation  was  between  7  and  02  with  a  value  of  -0.61. 

This  analysis  can  be  extended  by  looking  at  the  collinearity  of  all  the  parameters  used 
in  the  optimization  rather  than  simply  comparing  two  at  a  time.  The  collinearity  is  a  mea¬ 
sure  of  linear  dependencies  between  sets  of  parameters.  As  the  collinearity  value  increases 
the  parameters  are  more  related.  As  a  rule  of  thumb,  models  with  a  collinearity  value  less 
than  about  20  are  classified  as  identifiable  (Soetaert  and  Petzoldt  2010).  In  the  throm¬ 
bopoiesis  model,  the  collinearity  of  the  optimized  parameters  is  36.  In  the  granulopoiesis 
model,  the  collinearity  of  the  optimized  parameters  is  7.73,  and  in  the  lymphopoiesis  model, 
the  collinearity  of  the  optimized  parameters  is  7.80.  Thus,  the  optimized  values  for  the 
thrombopoiesis  model  should  be  interpreted  cautiously  because,  while  providing  accurate 
predictions  of  radiation  response,  they  may  not  provide  biologically  meaning  values. 
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Figure  C.l:  Thrombopoiesis  MCMC  analysis  results  showing  distributions  of 
optimized  parameters  and  correlations.  The  units  for  7  and  5  are  d-1,  and  the 
units  for  D®  are  Gy.  The  other  parameters,  A  and  02,  are  unitless.  The  parameter 
frequency  histograms  show  the  relative  likelihood  of  parameter  values  based  on 
how  often  they  occurred  in  the  MCMC  analysis. 
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Figure  C.2:  Granulopoiesis  MCMC  analysis  results  showing  distributions  of  op¬ 
timized  parameters  and  correlations.  The  units  for  ip  and  rj  are  d-1.  The  units 
for  I},,  D °00,  and  D?2  are  Gy.  The  other  parameter  l  is  unitless.  The  parameter 
frequency  histograms  show  the  relative  likelihood  of  parameter  values  based  on 
how  often  they  occurred  in  the  MCMC  analysis. 
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Figure  C.3:  Lymphopoiesis  MCMC  analysis:  Parameter  distributions  and  cor¬ 
relations.  The  units  for  7  are  d_1.  The  units  for  D®  are  Gy  and  02  is  unitless. 
The  parameter  frequency  histograms  show  the  relative  likelihood  of  parameter 
values  based  on  how  often  they  occurred  in  the  MCMC  analysis. 
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Abbreviations,  Acronyms,  and  Symbols 


ARA 

ARS 

0 

CFU-MK 

Co-60 

d 

DTRA 

FDA 

FIA 

FME 

G-CSF 

GM-CSF 

Gy 

GyEq 

7 

h 

HENRE 

HSCs 

IAEA 

IL-3 

IL-12 

kg 

MCMC 

METREPOL 

MKs 

MLT 

L 

n 

NETs 

NHP 

NK 

ODE 


Applied  Research  Associates,  Inc. 
acute  radiation  sickness 
beta  dose 

megakaryocyte  colony- forming  unit 

Cobalt  60 

day 

Defense  Threat  Reduction  Agency 
Food  and  Drug  Administration 
free- in- air 

flexible  modeling  environment  for  modeling,  sensitivity,  and  Monte  Carlo 
analysis  (R  plug-in) 

granulocyte  colony-stimulating  factor 
granulocyte  macrophage  colony-stimulating  factor 
gray 

gray  equivalent 
gamma  dose 
hour 

Health  Effects  from  Nuclear  and  Radiation  Environments 

hematopoietic  stem  cells 

International  Atomic  Energy  Agency 

interleukin  3 

interleukin  12 

kilogram 

Markov  Chain  Monte  Carlo 

Medical  Treatment  Protocols  for  Radiation  Accident  Victims  as  a  Basis  for 
a  Computerised  Guidance  System 

megakaryocytes 
midline  tissue 
microliter 
neutron  dose 

neutrophil  extracellular  traps 
non-human  primate 
natural  killer 

ordinary  differential  equation 


119 


R 

RBE 

RIPD 

Sv 

TPO 

UNSCEAR 

WMD 


software  programming  language  for  statistical  computing 
relative  biological  effectiveness 
Radiation-Induced  Performance  Decrement 
sievert 

thrombopoietin 

United  Nations  Scientific  Committee  on  the  Effects  of  Atomic  Radiation 
weapons  of  mass  destruction 
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